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A  CLOSED-FORM  DYNAMIC  ELASTICITY  SOLUTION  TO  THE  FLUID/ 
STRUCTURE  INTERACTION  PROBLEM  OF  A  TWO-LAYER  INFINITE 
VISCOELASTIC  CYLINDER  WITH  INNER  AND  OUTER  FLUID  LOADING 
SUBJECT  TO  FORCED  HARMONIC  EXCITATION 


INTRODUCTION 

This  research  treats  the  two-layer  cylinder  with  inner  and  outer  fluid  loading  subject  to  forced 
harmonic  vibration  within  the  scope  of  the  linearized  theory  of  dynamic  elasticity.  In  previous 
efforts,  investigations  of  the  statically  loaded  single-layer  cylinder  were  documented  by  numerous 
authors,  and  the  forced  vibration  response  of  the  single-layer  cylinder  was  investigated  by 
Holden,^  Kaul  and  McCoy, ^  and  others. 

In  the  theoretical  development  presented  in  this  report,  the  case  of  nonaxisymmetric  vibration 
is  discussed  in  terms  of  circumferential  order  number  n.  Results  are  presented  for  the  first  three- 
order  numbers,  namely,  n  =  0,  1,  and  2.  Three  excitation  states  are  considered:  the  normal 
pressure  P j,,  the  longitudinal  shear  stress  P^,  and  the  circumferential  shear  stress  Pq. 

The  development  proceeds  using  scalar  and  vector  displacement  potentials  to  describe  the 
motion  of  the  solid.  The  boundary  conditions  are  specified  in  terms  of  the  stresses  and 
displacements,  and  thus  must  be  rewritten  in  terms  of  the  displacement  potentials. 

For  the  case  of  the  two-layered  cylinder,  the  problem  is  reduced  to  solving  a  13-by-13  system 
matrix  for  the  undetermined  complex  coefficients.  This  matrix  requires  inversion  at  each 
location  in  the  wavenumber  and  frequency  plane,  allowing  all  dynamic  quantities  (stress,  strain, 
and  displacement)  to  be  calculated  from  the  coefficients.  Even  though  the  solution  is  closed- 
form,  a  sizeable  numerical  computation  is  still  necessary  to  extract  information  from  the 
governing  equations. 
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THE  PHYSICAL  MODEL 

The  physical  model  under  analysis  is  shown  in  figure  1.  A  fluid-filled  inner  cylinder 
(cylinder  1)  of  inner  radius  a  and  outer  radius  h  is  joined  to  a  second  cylinder  (cylinder 
2).  Cylinder  2  is  in  intimate  contact  with  cylinder  1  at  radius  r  =  h.  At  radial  distance  c, 
cylinder  2  is  in  contact  with  an  outer  fluid.  The  outer  fluid  extends  to  infinity  in  radial 
coordinate  r.  The  cylinders  and  fluids  extend  to  plus  and  minus  infinity  in  the  longitudinal  x- 
direction. 

As  stated  earlier,  the  system  is  treated  within  the  context  of  the  linearized  theory  of  elasticity. 
The  assumption  of  linearity  implies  that  both  the  cylinders  and  fluids  obey  linear  constitutive 
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equations.  For  the  cylinders,  this  means  that  there  is  a  linear  relationship  between  stress  and 
strain.  For  the  fluids,  deviations  from  the  hydrostatic  pressure  are  a  linear  function  of  the 
condensation  (fractional  density  change)  and  similarly  the  dilatation  (volume  strain).  Strains  are 
assumed  to  be  small  compared  to  unity  for  both  the  cylinders  and  the  fluids;  the  response 
therefore  is  linearly  related  to  the  excitation.  Both  media  are  restricted  to  small  signal  analysis, 
i.e.,  excitation  levels  on  the  acoustic  scale. 

Damping  is  included  by  the  use  of  a  complex  modulus  of  elasticity — the  well-known 
structural  damping  form.  The  fluids  are  modeled  as  inviscid  for  the  purpose  of  fluid  loading  the 
moving  structure.  The  radial  velocity  of  the  fluids  and  the  cylinders  is  identically  equal  atr  =  a 
and  r  =  c.  Under  the  inviscid  assumption,  the  fluids  interact  with  only  the  radial  motion  of  the 
cylinders  and  are  uncoupled  from  the  circumferential  and  longitudinal  motions.  A  perfect  slip 
condition  is  allowed  between  the  cylinders  and  the  fluids  for  the  longitudinal  and  circumferential 
directions.  This  assumption  allows  us  to  avoid  the  use  of  boundary  layer  physics  within  the 
context  of  this  development. 

As  described  in  the  Introduction,  three  harmonic  excitations  are  applied  by  the  outer  fluid  to 
cylinder  2  at  r  =  c:  radial  pressure  P„,  longitudinal  shear  stress  P,.,  and  circumferential  shear 
stress  Pq.  In  practice,  P^  and  Pq  would  only  be  developed  by  the  viscosity  of  a  real  outer 
fluid.  The  amplitude  distribution  (at  r  =  c)  of  the  excitations  in  the  wavenumber  and  frequency 
domain  must  be  known  a  priori',  otherwise  a  calculated  quantity  is  normalized  by  an  excitation 
stress  magnitude.  As  shown  in  figure  1,  the  excitations  are  at  a  definite  wavenumber  {k)  and 
frequency  (co)  and  are  of  circumferential  order  number  n.  The  response  of  the  cylinder  will 
therefore  be  a  function  of  k,  (O,  and  n,  as  well  as  of  the  geometric  and  material  properties  of  the 
composite  system. 
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EQUATIONS  OF  MOTION 
DISPLACEMENT  POTENTIALS 

The  equations  governing  the  motion  of  a  homogeneous,  isotropic,  linearly  elastic  body  are 
formulated  with  the  strain-displacement  relations  (equation  (1)),  the  constitutive  relations 
(equation  (2)),  and  the  stress  equations  of  motion  (equation  (3))  from  reference  3: 


^  =  (2) 

and 

(3) 

where  E,y  are  the  components  of  the  strain  tensor,  ty  are  the  components  of  the  stress  tensor,  m,  are 
the  components  of  the  displacement,  and  is  the  density  of  the  cylinder.  The  components  of 
acceleration  are  indicated  by  U-  and  the  body  forces  are  indicated  by  /b,.  The  commas  in 
equations  (1)  and  (3)  denote  partial  differentiation  as 


du- 

dxj 


dx- 


(4) 


where  the  components  of  the  vector  of  spatial  coordinates  Xj  are  equal  to  Xi  =  r,X2  =  0,  and  Xj  =  x. 

To  obtain  an  equation  only  in  terms  of  the  displacements  and  their  derivatives,  equations  (1) 
and  (2)  are  substituted  into  equation  (3),  resulting  in'* 


iK  +  \^c)\ji  =  PA- 


(5) 


Equation  (5)  is  expressed  in  vector  notation  as 


PA  «  +  +  l-^c)  VV  •  M  =  p  i/. 


(6) 
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The  displacement  vector  is  decomposed  with  the  scalar  and  vector  relationship 

M  =  V(p  +  Vx\|/  ,  (7) 

where  (p  =  (p  (r,  0,  x,  0  is  the  scalar  displacement  potential  and  \i/  =  V|/  (r,  6,  x,  t)  is  the  vector 
displacement  potential.  Equations  (6)  and  (7)  reduce  to 


For  the  present  cylindrical  geometry,  we  have  the  following  scalar  and  uncoupled  vector 


potentials: 

V2(p  =  \ip 

(12) 

and 

II 

(13) 
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Two  coupled  vector  potentials  in  r  and  0  are 


and 


r 


2^0 

,.230 


V72  ^0  2^¥, 


1 5  V, 

2-^  2 
c,dt 


1  ^  ¥e 

2..  2  ■ 
c,  dt 


(14) 


(15) 


The  Laplacian  is  defined  in  cylindrical  coordinates  as 


V 


2 


_3^  lA  ijl  a! 

'-dr  ^  rW^dx^' 


(16) 


Considering  the  scalar  potential  given  by  equation  (12)  and  the  harmonic  nature  of  the  forcing 
functions  in  figure  1,  we  will  assume  a  separation-of-variables-type  solution  in  r  and  0.  The 
scalar  potential  (p  (r,  0,  x,  t)  is  then  of  the  form 

tp  =  Q(r)0(0)e'^^"““'\  (17) 

Equation  (12)  yields  the  following  two  equations: 


and 


d  Q  1^ 
dr^  ^ 


f  0)^ 


2 

Q-^Q  =  0 


J  r 


+  nQ  =  Q. 


(18) 


(19) 


The  solution  of  equation  (18)  is  in  terms  of  Bessel  functions.  For  a  solid  cylinder,  the  field 
quantities  must  be  finite  at  the  center  (r  =  0);  therefore  only  Bessel  functions  of  the  first  kind, 
would  be  retained.  However,  for  hollow  cylinders,  such  as  cylinders  1  and  2,  Bessel  functions  of 
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both  the  first  and  second  kind  are  retained.  The  solution  to  equation  (19)  is  in  terms  of  the  sine 
and  cosine  of  argument  n9.  The  total  solution  is  therefore 


cp  =  [cos(n9)  +sin(n9)]  {A^t^ipr)  -<r  A2Y  ^{pr))  e  ^  ^  ^ 

with  p  defined  as 


(20) 


(21) 


A  similar  solution  for  the  A:-component  of  the  vector  potential,  is 

\|/^  =  [cos  (nd)  +  sin  (nQ)]  (B^J^iqr)  +  B2Y^{qr))  (22) 

with  q  defined  as 


2 

= 


(23) 


The  potentials  for  and  ij/e  can  be  uncoupled  by  noting  that  \|r^  and  \pe  contain  trigonometric 
functions  of  0  and  that  a  sine  dependence  on  9  in  \|/,.  is  equivalent  to  a  cosine  dependence  on  9  in 
\|/e  and  vice  versa.^  With  the  use  of  the  condition  V  •  \|/  =  0  to  eliminate  an  extra  constant,  \|/^ 
and  VI/9  become 


=  [Ci4+i(<?r)  +C2T„  +  i(^r)]  sin(nO)^'^^''  (24) 

and 

Ve  =  [-C/„+i(<?^) -C2T„^i(^r)]  cos(n0)e^'^^''  (25) 
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Displacement  Potential  Summary  for  Cylinders  1  and  2 

As  stated  earlier,  cylinders  1  and  2  require  Bessel  functions  of  the  first  and  second 
kind7  The  summary  of  the  scalar  and  vector  displacement  potentials  is 

^  +A2F^(pr)]  cos(n0)e'^^'*“®'\ 

+52F„(^r)]  sin(ne)e''^""~“'\ 

Ve  =  +-C2J'„  +  ,  (^r)]  cos(ne)e'^*''~®'\ 
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PARTICLE  DISPLACEMENT/DISPLACEMENT  POTENTIAL 


Intermediate  relationships  between  the  particle  displacements  and  the  displacement  potentials 
are  given  by  equations  (27)  through  (29),*  where  is  the  displacement  in  the  r-direction,  is  the 
displacement  in  the  circumferential  0-direction,  and  is  the  displacement  in  the  longitudinal  x- 
direction: 


^  r36  dx 


(27) 


and 


^  1^ 

dr 


_  ^ 

dx 


rdr 


1^¥,- 


(28) 


(29) 


STRESS/DISPLACEMENT  EQUATIONS 


In  the  problem  posed,  the  boundary  conditions  are  prescribed  with  the  stresses  and 
displacements.  Equations  (30)  through  (35)  relate  the  stress  components  to  the  particle 
displacements.  We  now  have  a  complete  progression  from  displacement  potentials  to  stresses, 
which  will  allow  us  to  solve  for  the  motion  of  the  composite  system: 


t 


rr 


dUc 

dr 


+ 


^  l^C- 
r  rdd  dx  ) 


+  2[ic 


d^c 

dr  ’ 


(30) 


'^ee 


=  X. 


du  u  i3v 

_ ^  4 - -4 - 

dr  r  rdQ 


dw^^ 
dx  j 


( 


+  2p,| 


r  rdd  I 


(31) 


X 


XX 


^  du,.  i3v 

,  _ ^  +  —  + - +  — 

‘■\dr  r  rdd  dx  ) 


+  2p. 


dWc 

dx 


(32) 
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T 


re 


^c\ 


dr 


r  rdd  ) 


(33) 


T 


ex 


( 


rdd 


+ 


dx 


(34) 


and 


^r.  = 


^  du  dw  '' 


(35) 


The  indicial  notation  is  interpreted  in  the  following  manner.  The  first  index  refers  to  the  plane 
normal  to  the  coordinate  identified  by  the  first  index.  The  second  index  refers  to  the  direction  in 
which  the  stress  acts.  Stresses  t,.,.,  Xqq,  and  are  normal  stresses,  with  x,.,.  acting  in  the  plane 
perdendicular  to  the  radial  coordinate  in  the  radial  direction.  The  remaining  three  stresses  (x,.0, 
Xqj.,  and  Xi-x)  are  shear  stresses,  with  x,.0  acting  in  the  plane  perpendicular  to  the  radial  coordinate  in 
the  circumferential  0-direction. 


STRAIN/DISPLACEMENT  EQUATIONS 

The  relationships  between  strain  and  displacement  in  radial  coordinates  are 


e 


rr 


dUc 
dr  ’ 


(36) 


^00 


l^c 

r  rdd 


(37) 


and 


(38) 
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STRESS,  STRAIN,  AND  DISPLACEMENTS  IN  TERMS  OF 
DISPLACEMENT  POTENTIAL  EQUATIONS 


The  stress  and  displacement  equations  needed  to  compose  a  suitable  set  for  the  application  of 
the  boundary  conditions  are  written  in  terms  of  the  displacement  potentials  by  substituting 
equations  (27)  through  (29)  into  equations  (30),  (33),  and  (35)  to  obtain  equations  (39),  (40),  and 
(41): 


=  A, 


(  2  2  2'' 
^  I  1^  9c  (Pc 

^  r dr  ■^r2ae2  '^dx^ 


+  2p 


d%c  l^V;,  ^Ve 


"1^3r2 


r^dQ  rdQdr  dxdr 


(39) 


!  2  2  2  2  2  ^ 

l^Vr  l^Ve  l^N^x  l^T 


r^dQ  ^  rdQdr  dxdr  rdx  rdxdQ  dr^ 


V 


^  rdr  '''r2302 


(40) 


X  =  LL 

rx 


2  2  2  2  2 

a  9c  1  laVe  l^ 

2.^-^  +  -—  — —  +—  -^¥e  +  -—  + 


V 


‘drdx  r^dQ  rdQdr  dx^  dr^ 


rdr  rdQdx 


(41) 


This  same  substitution  is  performed  on  strain  equations  (36),  (37),  and  (38)  to  obtain  equations 
(42),  (43),  and  (44): 


l^Vx  ^Ve 

dr^  r^dQ  rdQdr  dxdr  ’ 

2  2  2 
^  I^C  1^  9c  l^¥c  1^ 

^99  ^2d9  rdx  ^  r^dQ^  ^  rdxdQ  rdQdr  ’ 

2  2  2 

_  a  9c  ^  ¥e  1^¥9  la  ¥c 
dx^  dxdr  rdx  rdxdQ 


(42) 


(43) 


(44) 
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Equations  (45),  (46),  and  (47)  are  now  repeated  to  aid  the  reader: 


dr  r36 


(45) 


rd6  ^  dx  dr 


(46) 


and 


w 


C 


,  1  ,,  _l^r 

dx  dr  r30  ' 


(47) 
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OUTER  FLUID 

As  the  surface  of  the  cylinder  is  set  into  motion  by  the  excitation  of  a  pressure  wave  of 
magnitude  P^,  a  pressure  field  is  generated  in  the  outer  fluid  by  the  cylinder  surface.  When 
the  longitudinal  wavenumber  associated  with  the  vibration  of  the  cylinder  surface  is  smaller  than 
co/c^,  the  pressure  radiated  into  the  outer  fluid  is  in  the  form  of  a  propagating  wave.  If  the 
longitudinal  wavenumber  of  the  cylinder  vibration  is  greater  than  co/c^,  the  pressure  field  decays 
exponentially  from  the  surface  of  the  cylinder  in  the  positive  radial  direction.  The  outer  fluid  is 
modeled  as  an  ideal  (nonviscous)  linearly  elastic  media  that  cannot  sustain  shear  stress,  even 
when  it  is  in  motion.’  The  equation  of  state  is  then 

%  =  -Pfij.  (48) 


whose  components  are  given  by  Hooke’s  law  (a-j  =  )  by  setting  py  =  0: 


^11  “  ^22  “  ^33  “  ~  Pf> 

~  ^13  ~  ^23  “  ^  ■ 


(49) 


The  following  scalar  displacement  potential  describes  the  dynamic  motion  of  the  fluid  in 
cylindrical  coordinates:*® 


(Pf  = 


_  1  a  tp^ 


.29^2 


(50) 


where 


a^  la  1  a^ 


(51) 


In  the  case  of  fluids,  which  is  equivalent  to  where  p  is  the  bulk 

modulus  of  the  fluid.  The  particle  displacement  and  the  pressure  in  the  fluid  can  be  expressed  in 
terms  of  the  displacement  potential  as  follows: 

My  =  VtPy  (52) 
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d  (p, 

Pf  =  -P/T2  • 


For  propagating  wavenumbers  in  the  outer  fluid  field,  denoted  by  subscript  s,  the  scalar 
displacement  potential  is  given  by^^ 


aJ  iLr(l)  r  \  ]  /  i(kx-m) 

=  lcos(ne)c 


2  0)  ,  2  o 

^2  = 


00  ,2 
->k  . 


For  an  outgoing  wave  described  by  c  ’  ,  the  pressure  field  in  the  outer  fluid  is  physically  a 

radiation  from  the  surface  of  the  moving  cylinder.  The  radial  dependence  in  equation  (54)  is  a 
Hankel  function  of  the  first  kind,  indicated  by  superscript  (1);  i.e., 

«r(*2'')  =  (*2'-)  + (*2'-)  • 

Differentiating  equation  (54)  with  respect  to  radius  and  time  yields  equations  (55)  through 
(57)  for  the  fluid  pressure,  particle  displacement,  and  particle  acceleration,  respectively: 


jcos  (ne)e'^^"' 


s,  ^  .  a\  Hkx-(Ot) 

=  Jcos(n0)c 


2 

^  ^  r  rr(l)  /  \  ^  /  Cl\  i(.kx-(i>t) 

(^2'')  Jcos(ne)e 


For  nonpropagating  wavenumbers  in  the  outer  fluid,  the  evanescent  field  is  given  by  the  following 
displacement  potential,  where  the  radial  dependence  is  in  the  form  of  the  modified  Bessel 


function  K„: 


/X  N  /  iikx-m)  ,2  ,2  CO  for 

9^  =  HK^if^r)  cos  (nQ)  e  ,  =  k  -  —  , 

c. 


,2  CO 

k  >  — 
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Similar  to  the  above  discussion  of  the  propagating  case,  the  pressure,  particle  displacement,  and 
particle  acceleration  are  given  by  equations  (59)  through  (61)  as 


(f^r)  cos  (nd) 


(59) 


u 


s 


(f^r)  cos  (/t0)  e 


i  {kx-m) 


(60) 


and 


2 

8 


dfi 


COS(«e) 


(61) 


Figure  2  connects  the  wavenumber-frequency  plane  with  the  time-space  characteristics  of  the 
field  properties  in  the  outer  fluid.  The  vertical  axis  represents  wavenumber  with  respect  to  the 
longitudinal  x-axis  of  the  cylinder.  The  horizontal  axis  is  the  frequency  axis.  The  positive 
wavenumber  direction  corresponds  to  propagation  in  the  positive  x-direction,  and  the  negative 
wavenumber  direction  corresponds  to  propagation  in  the  negative  x-direction.  There  are  two 
distinct  regions  in  the  wavenumber-frequency  plane  corresponding  to  propagating  and  non¬ 
propagating  fields  in  the  outer  fluid. 

In  review,  variable  k  refers  to  the  wavenumber  that  exists  on  the  surface  of  the  cylinder.  The 
vibrating  cylinder  surface  will  initiate  a  propagating  field  in  the  outer  fluid  if  (sp-lc,^  >  k^.  When 
this  condition  exists,  the  radial  displacement  (m^)  of  the  cylinder  surface  must  project  (m^.’)  onto 
the  propagation  direction  of  the  field  in  the  outer  fluid.  For  a  given  frequency,  the  wavelength  of 
cylinder  vibration  must  be  greater  than  the  corresponding  wavelength  of  vibration  that  the  fluid 
will  propagate.  The  axisymmetric  n  =  0  vibration  produces  a  symmetric  field  in  the  outer  fluid 
as  a  function  of  0.  The  spatial  representation  of  a  point  (^i,C0i)  in  the  region  of  propagating 
wavenumber-frequency  pairs  given  by  |^^|>|^i|>0,  is  shown  in  figure  3,  where  Lj  is  the 
wavelength  and  is  the  velocity  of  propagation.  In  the  limit  of  |A:j|  ^  0 ,  the  radiated  field  is 
physically  a  series  of  cylindrical,  radially  spreading  wavefronts,  emanating  outward  from  and 
parallel  to  the  surface  of  the  cylinder.  The  intermediate  range  shown  in  figure  3  results  in  the 
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Figure  2.  Outer  Fluid  Propagating  and  Nonpropagating  Regions  of  the  Wavenumber- 

Frequency  Plane 
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wavefronts  propagating  out  along  conical  planes,  forming  an  angle  a  with  the  longitudinal  x-axis 
of  the  cylinder.  The  direction,  a,  of  the  radiated  field  will  be  given  by  equation  (62)  as 


a  =  acos 


^5  ^1 

acos—  =  acos—, 


(62) 


where 


,  2n  (0 

S  S 


(63) 


Propagation  ceases  when  This  condition  requires  that  c^  =  c^;  therefore,  the 

propagation  in  the  fluid  would  be  along  the  longitudinal  x-direction.  However,  there  is  only 
coupling  between  the  cylinder  motion  and  a  propagating  fluid  pressure  field  if  the  radial 
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displacement  u  can  be  projected  onto  the  propagation  direction  of  the  fluid  pressure  wave.  This 
is  not  possible  when  since  u  and  x  are  orthogonal  to  each  other;  this  condition  is  the 

onset  of  the  nonpropagating  region  shown  in  figure  2. 

The  Hankel  function  variation  of  the  outer  fluid  displacement  potential  can  be  observed  in  the 
inset  TIME-SPACE  diagrams  in  figure  2.  These  diagrams  depict  the  radial  spatial  variation  of 
the  magnitude  of  the  field  quantities  at  a  frozen  instant  of  time. 

As  a  line  of  constant  frequency  is  traversed  in  figure  2,  the  cylinder  wavelength  becomes 
shorter,  ultimately  reaching  a  length  less  than  the  length  of  the  corresponding  fluid  wavelength 
that  supports  propagation.  Under  this  condition,  there  is  no  longer  a  component  of  the  radial 
displacement  of  the  cylinder  that  can  be  projected  onto  the  propagation  direction  of  the  field  in  the 
outer  fluid;  equation  (62)  is  no  longer  valid;  and  the  field  in  the  outer  fluid,  p^,  induced  by  cylinder 
vibration  ceases  to  propagate.  The  radial  variation  of  the  field  quantities  is,  according  to  the 
modified  Bessel  function  depicted  in  figure  2,  for  the  nonpropagating  regions.  The  magnitude  of 
all  field  quantities,  such  as  the  displacement  potential  and  all  parameters  related  to  it,  undergo  the 
same  qualitative  spatial  variation. 
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INNER  FLUID 

Similar  to  the  behavior  of  the  outer  fluid,  as  the  surface  of  the  cylinder  is  set  into  motion  by 
the  excitation  of  a  pressure  wave  of  magnitude  P^,  a  pressure  field  (p,)  is  generated  in  the  inner 
fluid  by  the  cylinder  surface.  When  the  longitudinal  wavenumber  associated  with  the  vibration 
of  the  cylinder  surface  is  smaller  than  to/c„  the  pressure  radiated  into  the  inner  fluid  is  in  the  form 
of  a  propagating  wave.  If  the  longitudinal  wavenumber  of  the  cylinder  vibration  is  greater  than 
co/c/,  the  pressure  field  decays  exponentially  from  the  surface  of  the  cylinder  in  the  direction  of 
decreasing  radius. 

When  considering  the  behavior  of  the  field  at  r  =  0,  we  make  the  observation  that  and 
approach  infinity  as  the  argument  of  the  functions  approaches  zero.  However,  because  the 
interior  pressure  field  must  be  finite  at  r  =  0,  these  solutions  are  inappropriate  and  must  be 
discarded.  Bessel  functions  J „  and  /„  will  be  retained  to  describe  the  displacement  potentials  in 
the  inner  fluid. 


For  the  case  of  inner  fluid,  c.  =  where  p,-  is  the  bulk  modulus  of  the  fluid  and  p,  is 

the  density.  The  particle  displacement  and  the  pressure  in  the  fluid  can  be  expressed  in  terms  of 
the  displacement  potential  as  follows: 


and 


u^  =  Vtp. 


(64) 


2 

Pi  =  -Pi— 

dt 


(65) 


For  propagating  wavenumbers,  the  scalar  displacement  potential  is  given  by 


=  DJ^{g^r)co^{nd)e 


i(kx-(x)t) 


2  (0  ,2  f 

gi  =  ^-k  ,  for 


CO  ,2 

->k 

^i 


(66) 
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Differentiating  equation  (66)  with  respect  to  radius  and  time  yields  equations  (67)  through 
(69)  for  the  fluid  pressure,  particle  displacement,  and  particle  acceleration,  respectively: 


i{kx-  (Of) 


and 


p.  =  to^p.D/^  (gjr)  cos  (n0)  e 


M;  =  (g,r)  cos  (ne)  e 


^  T  r  \  /  iikx-m) 

—  =  ^  i8\r)  cos  (nQ)  e 


(67) 

(68) 


(69) 


For  nonpropagating  wavenumbers  in  the  inner  fluid  (imaginary  gi),  the  field  decays 
evanescently  and  is  described  by  the  following  displacement  potential,  where  the  radial 
dependence  is  in  the  form  of  the  modified  Bessel  function  /„,  evaluated  for  real  argument  m: 


(p.  =  Gl^  (mr)  cos  (n0)  e 


i  (Of) 


m 


,2  CD 
=  ^ 

C; 


for 


,2  © 

k 


(70) 


In  this  case,  /„  replaces  using  the  identity  /„  (x)  =  /  (ix) .  Similar  to  the  discussion  of 
the  propagating  case,  the  pressure,  particle  displacement,  and  particle  acceleration  are  given  by 
equations  (71)  through  (73)  as 


Pi 


Cd^p-GI^  (mr)  cos  (nQ)  e 


i  (Of) 


(71) 


G^I  (mr)  cos  (n0)  e 
dr  ” 


i  (kx~(i)t) 


(72) 


and 


2 

d  u- 


G^I„  (mr)  cos  (nQ)  e 
dr  " 


i  {kx-m) 


(73) 
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The  FORTRAN  algorithm  that  is  used  to  evaluate  the  Bessel  function  numerically  accepts 
complex  arguments.  Thus,  equations  (70)  through  (73),  which  use  Bessel  function  /„  to  replace 
as  the  argument  of  the  function  becomes  imaginary,  are  not  used  in  the  FORTRAN  algorithms 
developed  to  implement  the  solution  of  the  inner  fluid  displacement  potentials.  The 
development  presented  above  is  provided  in  the  event  that  the  complex  argument  Bessel  function 
routines  might  not  be  available  to  the  reader. 

Figure  4  connects  the  wavenumber-frequency  plane  with  the  time-space  characteristics  of  the 
field  properties  in  the  inner  fluid.  The  vertical  axis  represents  wavenumber  with  respect  to  the 
longitudinal  x-axis  of  the  cylinder.  The  horizontal  axis  is  the  frequency  axis.  The  positive 
wavenumber  direction  corresponds  to  propagation  in  the  positive  x-direction,  and  the  negative 
wavenumber  direction  corresponds  to  propagation  in  the  negative  x-direction.  There  are  two 
distinct  regions  in  the  wavenumber-frequency  plane  corresponding  to  propagating  and  non¬ 
propagating  fields  in  the  inner  fluid. 
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Figure  4.  Inner  Fluid  Propagating  and  Nonpropagating  Regions  of  the  Wavenumber- 

Frequency  Plane 


24 


TR  11,067 


SINGLE-LAYER  CYLINDER/FLUID  INTERFACES 
BOUNDARY  CONDITIONS 

The  solution  to  the  solid  cylinder  in  contact  with  an  outer  fluid  is  obtained  by  solving  for  the 
constants  Aj  ,^2  ,By  ,8^  ,Cj  and  in  equation  (26),  M  in  equation  (54),  and  D  in 
equation  (66)  for  the  pertinent  boundary  conditions  by  considering  the  state  of  stress  and 
displacement  at  the  cylinder  and  fluid  interfaces  r  =  a  ondr  =  b.  (Note  from  this  point  on  in 
the  derivations  that  subscript  c  (Lame  constants)  has  been  replaced  with  1  or  2  to  indicate  cylinder 
1  or  2.  Superscript  Cl  or  C2  has  been  added  to  all  other  pertinent  quantities  to  distinguish  the 
first  from  the  second  cylinder.) 

The  final  problem  contains  a  total  of  seven  unknown  constants.  There  are  seven  boundary 
conditions,  which  are  reflected  in  rows  1  through  7,  respectively,  of  system  matrix 
smc\.  Initially  the  problem  was  set  up  with  eight  boundary  conditions,  where  the  eighth  equated 
the  radial  displacement  of  the  cylinder  to  the  outer  fluid  dX  r  =  b.  However,  due  to  the 
properties  of  the  cylinder  used  for  the  trial  problem  and  the  limited  precision  available  for  the 
calculation,  rows  4  and  8  of  this  original  matrix  were  not  sufficiently  different  to  enable  the 
required  inversion,  a  condition  refered  to  as  algorithmically  singular. 

To  circumvent  this  original  inversion  problem,  the  radial  displacement  boundary  condition 
evaluated  at  r  =  b, 


Cl 


r  =  b 


dr 


(74) 


was  assembled  and  written  in  terms  of  the  unknown  fluid  constant  M,  assuming  co^/c^  >  ,  as 

follows: 


M  = 


r(l) 


+  Cj  ikJ +  C2  ikY^^^iq-^b)  I- 


(75) 
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Equation  (75)  was  used  in  the  assembly  of  boundary  condition  1  (equation  (76))  to  eliminate  the 
outer  fluid  constant  M.  This  added  a  fluid  term  to  each  of  the  terms  in  boundary  condition  1, 
making  it  sufficiently  different  from  boundary  condition  4  to  allow  the  numerical  inverse  of 
matrix  smc\  to  exist.  Thus,  the  system  was  reduced  from  eight  equations  for  the  boundary 
conditions  to  seven,  which  also  decreased  the  time  needed  for  computation. 

The  seven  boundary  conditions  are  given  by  equations  (76)  through  (82)  as 


Boundary  condition  1 


Boundary  condition  2 


Boundary  condition  3 


Boundary  condition  4 


Boundary  condition  5 


Boundary  condition  6 


Boundary  condition  7 


Cl 


r  =  b 


(-I’c-Ps) 


Cl 


r  =  b 


Cl 

'fre 


r  =  b 


Cl 


(76) 

(77) 

(78) 

(79) 

(80) 

(81) 

(82) 
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The  fourth  boundary  condition  (equation  (79))  for  radial  stress  is  assembled  by  substituting 
equation  (67)  and  the  appropriate  derivatives  of  equation  (26)  into  equation  (39).  The  resulting 
expression  is  in  terms  of  the  unknown  constants  for  the  displacement  potentials  and  is  given  as 

Cir2|a,n/  3  >,■]  3  V 

+  +  )cos  (nG)  e' 

=  - a)2p.D/„  (g^a)  cos  (nG)  .  (83) 
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As  previously  described,  the  first  boundary  condition,  equation  (76),  is  assembled  by  using 
equation  (75)  to  eliminate  the  outer  fluid  constant  M,  resulting  in 


(S2I’) 


2 

U 


2  f  \  f  2 

+A^‘  (\  +  2n,)^y„0,i>)+l-y>,i)  ,  -x,y>,4^  +  «:" 

a;b^.  (Sa*) 


2Pi»('  3 


b  d  [  „(i)  ,  T 
^  (^2^) 


Cl  2p,n/  3  \ 


"py  /fl^dib) 


Cl  .,  a 


2  «7’(«2'>) 


(«2^) 


i4+|(9l*) 


+ cf  2p,rt|-y„ , ^  ^y„ ,,(«,*) 

a;  (*2*) 


COS  (rt9)  =0.  (84) 
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SYSTEM  MATRIX 

The  seven  boundary  conditions  in  equations  (76)  through  (82)  are  assembled  in  matrix  form  in 
equation  (85): 


5mcl(l,l)  smc\{l,2)  smc\{\,'i)  smc\{\,4)  5mcl(l,5)  5mcl(l,6)  0 

'-Po 

smc\(2,\)  smc\{2,2)  smc]  (2,3)  smc\{2,A)  smc\  {2,5)  5mcl(2,  6)  0 

a‘^' 

-p^ 

5mcl  (3, 1)  ^mcl(3,  2)  5»icl(3,  3)  smc\  (3,4)  jwcl(3,  5)  jmcl(3,  6)  0 

-Pb 

^mcl  (4,  1)  smc\{4,2)  smc\ {4,3)  smc\{4,4)  5mcl(4,  5)  smc\{4,6)  smc\{\,l) 

= 

0 

$mc\{5,\)  smc\{5,2)  5wcl(5,  3)  imcl(5,  4)  5mcl(5,  5)  5mcl  (5,  6)  0 

cr 

0 

5wcl  (6,  1)  smcl{6,2)  smcl{6,3)  ^mcl(6, 4)  5mcl(6,  5)  ^mcl(6,  6)  0 

0 

smc\{7,\)  smc\{l,2)  smc\{l,3)  smc\{l,4)  smc\ {7,5)  smc\ {7,6)  smc\{],7)^ 

D 

0 

(85) 


Elements  of  the  smc\  matrix  are  shown  in  equations  (86)  through  (134).  Solution  of  equation 
(85)  at  each  point  in  the  wavenumber-frequency  plane  for  the  vector  of  unknown  constants  allows 
the  displacement  potentials  of  the  system  to  be  evaluated  with  equations  (26),  (54),  and 
(58).  Equation  (85)  will  be  evaluated  for  only  one  nonzero  forcing  function  at  a  time.  The 
coefficients  will  be  normalized  by  the  forcing  function  stress,  which  will  cast  the  output  in  terms 
of  a  transfer  surface  for  the  particular  cylinder  variable  chosen,  i.e.,  displacement,  velocity,  stress, 
or  strain. 


From  boundary  condition  1, 


5mcl(l,  1)  =  (?^i +  2Pi)^/„(Pih)  +  ^/„(Pih) 


_J  +  p  CO  - - 

b  d  (  „(l) 


V 


(Sob) 


(86) 
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X 


(^2^) 


-\yrfPxb)[l^^k^ 


(87) 


,smcl(l,3)  = 


2  (1) 

<*2*)  j 


(88) 


smcl  (1,4)  =  -: 


(89) 


^mcl  (1,  5)  = 


^(1)  . 

l(9,(>)  +  (<:Py  g/  "  ^ |(9l*) . 

3?l""  (*2*> 


(90) 


smcl  (1,  6)  =  i2p,c"-K„,  ,(g,(,)  +  iCpy-g,  "  ^ - ,^„  ,(?,*) , 

57l"«  <«2«>) 


(91) 


and 


smcl  (1,  7)  =  0. 


(92) 


From  boundary  condition  2, 


5mcl(2,  1)  =  i2\i^k^J^{p^b),  (93) 

mcl(2,  2)  =  i2\x.^k^YJ^^b),  (94) 
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i\)i,kn 

smc\{l,2))  = 


smc\  (2, 4) 

smc\  (2,  5) 

smcl  (2,  6) 


(95) 


(96) 


and 


smcl  (2,  7)  =  0. 


(99) 


From  boundary  condition  3, 


smc\  (3,  1) 

= 

(100) 

smcl  (3,  2) 

=  1' 

(101) 

5mcl  (3,  3) 

(102) 

smcl  (3,  4) 

(103) 

smcl  (3,  5) 

(104) 

smcl  (3,  6) 

(105) 

and 
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smcl  (3,7)  =  0.  (106) 

From  boundary  condition  4, 

smcl  (4,1)  =  +  (107) 

.mcl(4,2)  =  (Xj  +  2iij)^F„(Pja)  +  ^|.F„(Pja)-^V„(p,a)-^2^,F„(p,a),  (108) 


smcl  (4,  3)  =  -2|X,^(^/„(^ja)  -a|-/^(<7ja)j,  (109) 

smcl{4,4)  =  -^l^i^[Y^(q^a)-a~Y„(q^a)y  (HO) 

5mcl(4,5)  =  2^^ik^J^^^(q^a),  (HI) 

^mcl(4,6)  =  2[i^ik^Y^_^^(q^a),  (112) 

and 

^mcl(4,7)  =  G)2p.£>/^(^ja).  (H3) 


From  boundary  condition  5, 

smcl(5,l)  =  i2\L^k^J^(p^a),  (114) 

smcl  (5,2)  =  /2|aj^|-F„(pja),  (115) 

ili,kn 

smcl(5,3)  =  ——J^iq^a),  (116) 
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i\y.,kn 

smc\{5,A)  =  -^YJ,q^a),  (117) 

a  "  ‘ 

smcia.i)  =  (/,.,(«, (118) 

jmcl(5,6)  =  H,(i'„,i(<?,<!)(^-<:^)-^l',+  l(?i<i)(^)-j^l'„  +  l('?lO)).  (IW) 
and 

5mcl(5,7)=0.  (120) 


From  boundary  condition  6, 

2|X,n/i  ^  \ 

smc\{6,\)  =—[-/>!«) -^/„(Pia)J,  (121) 

™cl(6,2)  =  (122) 

smcl(6,3)  =  +  (123) 

smc\(6,A)  =H,(-g^l',(9,'!)  +  ^^l’„(?l‘')-^l'„(«lO)).  (124) 

s»icl(6,5)  =  «H,(^/„+,(«,a)-/„»l(9l<')(^)).  (125) 

smc\(6,6)  =  «:l>i(|;l'„.,(9,<i)-l',+  ,(?,«)(i^)).  (12® 
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and 


smc\  (6, 7)  =  0. 

(127) 

From  boundary  condition  7, 

0 

^mcl(7,  1)  =^7„(Pja), 

(128) 

smc\{l,2)  =^YJ^^a), 

(129) 

smc\  (7,  3)  =  -J  (q^a), 
a  "  ^ 

(130) 

smcia,4)  =  -Y^iq^a), 
a  "  ^ 

(131) 

smcl  (7,5)  =  ikJ^_^^(q^a), 

(132) 

smcl  (7,6)  =  ikY^^^iq^a), 

(133) 

and 

,smcl(7,7)  =  ^^/„(gia). 

(134) 
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TWO-LAYER  CYLINDER/FLUID  INTERFACES 
BOUNDARY  CONDITIONS 

In  this  section,  we  will  formulate  the  13  boundary  conditions  required  for  the  two-cylinder 
problem  shown  in  figure  1 .  At  radius  r  =  c  (the  second  cylinder/fluid  interface),  the  radial  stress 
in  cylinder  2  is  set  equal  to  the  negative  of  the  magnitude  of  the  pressure  in  the  fluid.  Also  at 
radius  r  =  c,  the  radial  displacement  of  the  cylinder  equals  the  radial  displacement  of  the 
fluid.  At  radius  r  =  b,  there  is  a  continuity  of  displacements  and  stress  components  between 
each  cylinder  in  the  plane  of  contact  perpendicular  to  the  radial  coordinate.  At  radius  r  =  a  (the 
first  cylinder/fluid  interface),  the  radial  stress  in  cylinder  1  is  set  equal  to  the  negative  of  the 
pressure  in  the  inner  fluid.  Additionally,  the  radial  particle  displacement  of  the  fluid  and 
cylinder  is  equated.  The  shear  stresses  and  vanish  due  to  the  inviscid  assumption.  For 
the  composite  problem,  the  following  description  at  the  boundaries  exists: 


Boundary  condition  1 

C2 

^rr 

=  ^-Pn-P.)\ 

(135) 

\r  ~  c 

\r  =  c 

Boundary  condition  2 

Cl 

=  -Px 

, 

(136) 

r  =  c 

\r  =  c 

Boundary  condition  3 

Cl 

=  -^0 

, 

(137) 

r  =  c  ^ 

\r  =  c 

Boundary  condition  4 

Cl 

'^rr 

Cl 

(138) 

r  =  ^  ^ 

r  =  h 

Boundary  condition  5 

Cl 

^rx 

Cl 

(139) 

r  =  h 

r  ~  h 

Boundary  condition  6 

Cl 

Cl 

(140) 

r  =  b 

r  =  h 

Boundary  condition  7 

Cl 

Cl 

r  =  h  ^ 

r  =  b 

(141) 

Boundary  condition  8 

Cl 

II 

il 

< 

r  =  b 

(142) 
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Boundary  condition  9 

C2 

II 

11 

r  =  b 

(143) 

Boundary  condition  10 

Cl 

=  -Pi\ 

r  -  a  Ir 

=  a 

(144) 

Boundary  condition  11 

Cl 

X 

rx 

=  0. 

r  =  a 

(145) 

Boundary  condition  12 

Cl 

^.0 

=  0. 

r  =  a 

(146) 

Boundary  condition  13 

Cl 

=  “/I 

r  =  a  Ir  = 

-  a 

(147) 

The  same  procedure  employed  for  the  single-layer  cylinder,  which  combined  the  outer  fluid 
radial  displacement  condition  with  the  boundary  condition  for  radial  stress,  is  used  in  the  two- 
layer  cylinder  problem  formulation.  At  radius  r  —  c,  the  continuity  of  radial  displacement 
between  the  cylinder  and  the  outer  fluid  requires  that 


Cl 


u 


C 


(148) 


This  boundary  condition  is  assembled  with  equations  (26),  (27),  and  (54).  Solving  for  M  as 
follows  yields 


M  = 


r(l) 


(149) 


Equation  (149)  is  combined  with  equation  (135)  to  eliminate  the  constant  M,  reducing  the  size  of 
the  system  matrix  from  14  by  14  to  13  by  13. 
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The  13  boundary  conditions  correspond  to  rows  1  through  13,  respectively,  of  equation 
(150).  The  elements  of  the  system  matrix,  C,  are  given  in  equations  (151)  through 
(319).  Solution  of  equation  (150)  at  each  point  in  the  wavenumber-frequency  plane  for  the 
vector  of  unknown  constants  allows  the  displacement  potentials  of  the  system  to  be  evaluated 
with  equations  (26),  (54),  and  (58).  Matrix  C  will  be  evaluated  for  only  one  nonzero  forcing 
function  at  a  time  in  any  given  simulation.  The  coefficients  will  be  normalized  by  the  forcing 
function  stress  magnitude  (Pf^,  P^,  or  Pq),  which  will  cast  the  output  in  terms  of  a  transfer  surface 
for  the  particular  cylinder  or  fluid  variable  chosen,  i.e.,  displacement,  velocity,  stress,  strain,  or 
fluid  pressure. 

The  13  boundary  conditions  are  shown  diagrammatically  in  figure  6  for  the  two-layer  solid 
cylinder  immersed  in  an  outer  fluid.  All  boundary  conditions  are  in  effect  over  the  entire  surface 
of  contact  for  the  media  at  radii  r  =  a,b,  and  c. 
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Figure  6.  Boundary  Conditions  for  Two-Layer  Infinite  Cylinder  and  Fluids 
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SYSTEM  MATRIX 

As  noted  earlier,  the  13  boundary  conditions  are  assembled  into  matrix  form  in  the  following 
expression: 


C(l,l)  C(l,2)  C(l,3)  C(l,4)  C(l,5)  C(l,6) 
C(2,  1)  C(2,2)  C(2,3)  C(2,4)  C(2,5)  C(2,6) 
C(3,  1)  C(3,2)  C(3,3)  C(3,4)  C(3,5)  C(3,6) 
C(4,  1)  C(4,2)  C(4,3)  C(4,4)  C(4,5)  C(4,6) 
C(5, 1)  C(5,2)  C(5,3)  C(5,4)  C(5,5)  C(5,6) 
C(6, 1)  C(6,2)  C(6,3)  C(6,4)  C(6,5)  C(6,6) 
C(7,  1)  C(7,2)  C(7,3)  C(7.4)  C(7,5)  C(7,6) 
C(8,  1)  C(8,2)  C(8,3)  C(8,4)  C(8,5)  C(8,6) 
C(9,  1)C(9,2)  0  0  C(9,5)  C(9,6) 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 


0  0  0  0 

0  0  0  0 

0  0  0  0 

C(4,7)  C(4,8)  C(4,9)  C(4,  10) 

C(5,7)  C(5,8)  C(5,9)  C(5,  10) 

C(6,7)  C(6,8)  C(6,9)  C(6,  10) 

C(7,7)  C(7,8)  C(7,9)  C(7,  10) 

C(8,7)  C(8,8)  C(8,9)  C(8,  10) 

C(9,7)  C(9,8)  0  0 

C(I0,7)  C(10,8)  C(10,9)  C(10,  10) 
C(ll,7)  C(ll,8)  C(ll,9)  C(ll,  10) 
C(12,7)  C(12,8)  C(12,9)  C(12,  10) 
C(13,7)  C(I3,8)  C(13,9)  C(13,  10) 


0  0  0 

0  0  0 

0  0  0 

C(4,  II)  C(4,  12)  0 

C(5,  11)  C(5,  12)  0 

C(6,  11)  C(6,  12)  0 

C(7,  11)  C(7,  12)  0 

C(8,  11)  C(8,  12)  0 

C(9, 11)  C(9,  12)  0 

C(10,  II)  C(10,  12)  C(10,  13) 
C(ll,  11)  C(ll,  12)  0 

C(12,  11)  C(12,  12)  0 

C(I3,  11)  C(13,  12)  C(13,  13) 


-Po 

-P^ 

0 

‘-I 

0 

0 

= 

0 

a‘^' 

0 

0 

4' 

0 

c- 

0 

c^' 

0 

D 

0 

150) 


The  elements  of  equation  (150)  are  given  by  equations  (151)  through  (319).  From  boundary 
condition  1,  the  first  row  yields 


C(l,  1)  =  + 


2 

7  +  P,“ 


C(l,2)  = 


2 


'h  2  ^ 

-  +  p,co 


V 


dr\  ^ 


(151) 


(152) 
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C(l,3) 


C(l,4) 


2 

np^co 
^  c  3 


ar 


H 


(1) 


H 


(1) 


(^2^) 

(^2^) 


+ 


«p  (0 


Hl"  («,c) 


dr 


«y’  (S2^) 


;/„(^2^)’ 


Y^iq^c). 


C(l,5)  =  +  ^  Vn+l^‘^2^^’ 

(^2^) 


C  ( 1,  6)  =  ^  j(^2^)  +  a 


2 


^n  +  l('?2^)’ 


0(1,7)  =  0, 

0(1, S)  =  0, 

0(1,9)  =  0, 

0(1,  10)  =  0, 

0(1,  11)  =  0, 

0(1, 12)  =  0, 

[ 

0(1,  13)  =  0. 


(153) 

(154) 

(155) 

(156) 

(157) 

(158) 

(159) 

(160) 

(161) 

(162) 

(163) 
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From  boundary  condition  2, 


C(2,  1)  =  (164) 

C(2,2)  =  (165) 

i\i.ykn 

C(2,  3)  =  (166) 

iiirfkn 

C(2,4)  =  -^y,(?2c).  (167) 


C(2,  5)  ii2(^+l(^2^)(  ^2  (168) 

C(2,6)  =  (169) 


C(2,7)  =  0,  (170) 

C(2,8)=0,  (171) 

C(2,9)  =  0,  (172) 

C(2,10)=0,  (173) 

C(2,  11)  =0,  (174) 

C(2,  12)  =  0,  (175) 

ind 


C(2,  13)  =  0. 


(176) 
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From  boundary  condition  3, 


(177) 

"  1  [fn(P2^'> 

(178) 

C(3,3)  =  + 

(179) 

C(3,4)  =  [^2[-^Y^(q^c)  +  -^^Y^iq^c)-'^Y^iq^c)]^ 

(180) 

C(3,5)  -  «H2(|;7„,(?2c)-7„,(?2c)(  ^  j). 

(181) 

C(3,6)  -  z/:it2(^l^„  +  i(^2^)-l^n+i(^2^)(  c  )]> 

(182) 

C(3,7)  =  0, 

(183) 

C(3,8)  =  0, 

(184) 

C(3,9)  =  0, 

(185) 

C(3,  10)  =  0, 

(186) 

C(3,  11)  =  0, 

(187) 

C(3,  12)  =  0, 

(188) 

and 

C(3,  13)  =  0. 

(189) 
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From  boundary  condition  4, 

C(4,1)  =  +  + 

C(4,2)  =  +  + 

C(4,3)  = 

C(4,4)  =  -2n,^[YJg,b}-b^YJ(,,b)]. 

<2(4,5)  =  i2\Ljc^J^^^iq^b), 

C(4,6)  = 

C<4,7)  =  -a,+2n,)^/>,fr)-^l./.(p,i,)  + V„(Pl(’)(p  +  '^0' 

<2(4.8)  =  -(^  +  2n,)j^l'>,i)-^|:y>,6)  +  ^,r>,(,)[^  +  ^2j, 

C(4,9)  =  2n,A(7„(«,(,)-;,|-7.(,,i,)), 

C(4,10)  =  2n,^(y„(«,6)-i|-y,(,,i,)], 

C(4,  11)  =  -i2n,k^J^,,(q,b), 

C(4,  12)  = -;2n,l:|-F„^, (,,(.), 
and 

C(4,  13)  =  0. 


(190) 

(191) 

(192) 

(193) 

(194) 

(195) 

(196) 

(197) 

(198) 

(199) 

(200) 

(201) 

(202) 
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From  boundary  condition  5, 


C(5,l)  =  (203) 

C(5,2)  =  (204) 

i\iLr,kn 

C(5,3)  =  -|-/„(^2^),  (205) 

iu,^kn 

C(5,4)  =  -^Y^iq^b),  (206) 


C(5,5)  (207) 

C(5,6)  =  (208) 


C  (5,  7)  =  -/2|Xi^^4(Pjb) ,  (209) 

C(5,8)  =  -ilii^k^Y^ip^b),  (210) 

/Li,^n 

C(5,9)  = - (211) 

iix.kn 

C(5,  10)  = - '^Y^iq^b),  (212) 


C(5,ll)  (,,),)],  (213) 

/  2  \ 

C(5,12)  =-n,(r^^,(^,fc)(2^-t2j-|-K„^,(«,),)[!i|i)-AK^^l(,,i,)J,  (214) 

and 

C(5,  13)  =0.  (215) 
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From  boundary  condition  6, 

C(6,l)  =  (216) 

2)‘i''/'  1  fl  'l 

C(6,2)  =  (217) 

=H2(-3^-f„('?2*)  +  j|:2„(<72«-^/,(?2fc)),  (218) 

c<6,4)  +  am 

12(6,5)  =  |(^26)“2„,  ,(^2(')(^  ^  (220) 

C(6,6)  (221) 

C(6,7)  =-- ^(^-/„0,5)-|/>,6)J,  (222) 

C(6,8)  =  -^(^ly^(p,b)-±y^(p,l,)].  (223) 

(2(6,9)  =  “)(i(“^^2„(^j6)  +  ^^/,,(,7j5)- — 7„(^j(>)j,  (224) 

C(6,  10)  =  -^,(-^yjf^l,)^l^y^(^^l,)~^yj^m)-  (225) 

C(6,  11)  =  -«H,(|:/.+  ,(«,«-/„,(<(,6)(i^)),  (226) 

C(6,  12)  =  -ttn,(|-F,,,(,,i)-y„^,(,S,,6)[l±i!]),  (227) 

and 

C(6,  13)  =  0.  (228) 
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From  boundary  condition  7, 


CO  A) 

(229) 

COO)  =^yMb). 

(230) 

COO)  = 

(231) 

CO  A)  = 

(232) 

C{1,5)  =  + 

(233) 

C(7,6)  =  ikY^^^iq^b), 

(234) 

C(7,7)  =  -|-(/„(pib)), 

(235) 

C(7,8)  =-^Y^(p,b)), 

(236) 

C(7,9)  =  -lj„iq,b). 

(237) 

C(7,10)  = 

(238) 

C(7,  11)  = -ikJ^^^iq^b), 

(239) 

C(7,  12)  =  -/^F„  +  i(^ib), 

(240) 

and 

C(7, 13)  =0. 

(241) 
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From  boundary  condition  8, 


C(8,l) 

(242) 

C(8,  2) 

(243) 

C(8,3) 

(244) 

C(8, 4) 

(245) 

C(8,5) 

(246) 

C(8,  6) 

(247) 

C(8,7) 

II 

(248) 

C(8,8) 

(249) 

C(8,  9) 

(250) 

C(8,  10) 

(251) 

C(8,  11) 

(252) 

C(8,  12) 

(253) 

and 

C(8,  13)=  0.  (254) 
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From  boundary  condition  9, 


C(9, 1)  = 

(255) 

C(9,2)  =  ikY^ip^b), 

(256) 

C(9,3)  =  0, 

(257) 

C(9,4)  =  0, 

(258) 

C(9,5) 

(259) 

C(9,6) 

(260) 

C(9,7)  =  -ikJ^ip^b), 

(261) 

C(9,8)  = 

(262) 

C(9,9)  =  0, 

(263) 

C(9, 10)  =  0, 

(264) 

(265) 

C(9,12)  =|:i'„*i(9,6)  +  "j'i',,,(,,6). 

(266) 

and 

C(9, 13)  =0. 

(267) 
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From  boundary  condition  10, 


C(10,l)  =  0, 

(268) 

C(10,2)  =  0, 

(269) 

C(10,3)  =  0, 

(270) 

C(10,4)  =  0, 

(271) 

C(10,5)  =  0, 

(272) 

C(10,6)  =  0, 

(273) 

C(10,7)  =  (>.j  +  2n,)A7^(p,a)  +  -j|:J>ia)-V,(P,‘')( 

/ 

(274) 

2  ^ 

C(10,8)  =  {\  +  2\i^)  ^^^Y^{p^a)  +  -^Y^{p^a)-X^Y^{p^d)\ 

(275) 

C(10,9)  =  -l\i,^lj^{q,a)-ay^{q,a)). 

(276) 

C(10,10)  = 

{111) 

C(10,  11)  = 

(278) 

C(10,12)  =  i2\i^k^Y^^^{q^a), 

(279) 

and 

C(10,13)  =  p.co^D/„(gja). 

(280) 
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From  boundary  condition  11, 

C(ll,  1)  =  0,  (281) 

C(ll,2)  =  0,  (282) 

C(ll,3)  =  0,  (283) 

C(ll,4)  =  0,  (284) 

C(ll,5)=  0,  (285) 

C(ll,6)  =  0,  (286) 

C{n,l)  =  (287) 

C(ll,8)  =  (288) 

C{\\,9)  (289) 

(290) 


2  N 

C(ll.ll)  =  H, (/„,(<?, (291) 

2  \ 

C(ll,12)  =  +  (2M) 

and 

C(ll,  13)=  0.  (293) 
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From  boundary  condition  12, 

C(12,l)  =  0, 

C(12,2)  =  0, 

C(12,3)  =  0, 

C(12,4)  =  0, 

C(12,5)  =  0, 

C(12,6)  =  0, 

C(12.8)  = 

y  2 

C(12,9)  =  + 

C(12,10)  =  + 

C(12,U)  = 

C(12,12)  = 
ind 

C(12,  13)  =  0. 


(294) 

(295) 

(296) 

(297) 

(298) 

(299) 

(300) 

(301) 

(302) 

(303) 

(304) 

(305) 

(306) 
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From  boundary  condition  13, 

C(13,  1)  =0,  (307) 

C(13,2)  =  0,  (308) 

C(13,  3)  =  0,  (309) 

C(13,4)  =0,  (310) 

C(13,5)  =  0,  (311) 

C(13,6)  =0,  (312) 

c (13.  7)  a),  (313) 

C ( 13,  8)  =|:y„(p, a).  (314) 

C(13,9)  =  ^/„(^ia),  (315) 

C(13,10)  =  ^F„(^i«),  (316) 

C(13,ll)  =  (317) 

C(13,12)  =  ikY^^^{q,a),  (318) 

and 

C(13,13)  (319) 
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SOLUTION  METHODOLOGY 

As  discussed  in  the  previous  sections,  the  system  matrices  in  equations  (85)  and  (150)  provide 
the  undetermined  constants  that  allow  for  the  evaluation  of  the  displacement  potentials  directly 
from  equations  (26),  (54),  (58),  (66),  and  (70).  At  this  point  we  can  evaluate  all  of  the  dynamic 
quantities  of  the  composite  systems,  i.e.,  the  single-  or  two-layer  cylinders  with  fluids.  At  any 
cylinder  radius,  the  stresses  can  be  evaluated  with  equations  (39)  through  (41),  the  strains  with 
equations  (42)  through  (44),  and  the  displacements  with  equations  (45)  through 
(47).  Additionally,  the  inner  fluid  pressure  field  can  be  evaluated  with  equations  (67)  and  (71), 
and  the  outer  fluid  pressure  field  with  equations  (55)  and  (59).  Fluid  velocities  can  be  evaluated 
in  a  similar  fashion.  The  Bessel  functions  are  evaluated  with  the  series  representations  given  in 
appendix  A,  and  the  necessary  derivatives  are  formed  with  the  relations  given  in  appendix 
B.  The  equations  used  for  calculating  the  output  quantities  are  given  in  appendix  C.  The 
output  quantity  is  normalized  by  the  particular  excitation  chosen  (equations  (85)  and  (150)  are 
evaluated  for  only  one  nonzero  excitation  at  a  time).  Normalizing  the  output  quantity  by  the 
excitation  casts  the  results  in  a  transfer  surface  when  computation  is  over  a  range  of  wavenumbers 
and  frequencies. 

The  two-layer  model  is  compared  against  the  single-layer  model  for  an  equivalent  cylinder  in 
appendix  D.  The  two  simulations,  which  are  shown  overlaying  each  other,  provide  a  verification 
of  the  agreement  between  the  elasticity  models. 


NUMERICAL  CONSIDERATIONS 

The  procedure  described  above  is  repeated  over  the  domain  of  wavenumber  and  frequency  at 
regular  intervals.  The  interval  used  for  the  simulations  extending  to  5000  or  8000  Hz  is  20.0  Hz 
and  the  interval  used  for  the  60,000-Hz  simulations  is  200  Hz.  The  wavenumber  interval  used 
for  the  5000-Hz  simulations  is  0.257  rad/m;  1.545  rad/m  is  used  for  all  other  simulations.  These 
intervals  produce  transfer  surfaces  ranging  from  32,000  to  194,000  points  each. 

When  sharp  resonances  are  present,  a  visual/numerical  artifact  is  apparent  in  the  simulations 
due  to  the  finite  wavenumber  and  frequency  interval  used  for  evaluating  the  closed-form 
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solution.  This  is  most  evident  in  the  wire  frame  surface  representation  of  the  simulations,  which 
is  shown  later  in  the  report.  Small  asymmetries  with  respect  to  wavenumber  are  evident  in  the 
magnitude  cuts  as  well.  The  first  problem,  which  arises  when  a  Cartesian  wire  frame  surface 
algorithm  is  used  with  sloped  topography  in  the  xy  plane,  is  inevitable  because  of  the  graphics 
software  that  is  used  to  display  the  wire  frame  surfaces.  The  second  problem  results  when,  for  a 
given  frequency,  slightly  different  wavenumber  locations  are  being  calculated  in  the  minus  and 
plus  wavenumber  half-planes.  This  problem  could  be  minimized  by  decreasing  both  the 
frequency  and  wavenumber  interval  by  several  orders  of  magnitude,  although  a  corresponding 
increase  in  calculation  time  would  result.  The  image  sizes  were  selected  as  an  appropriate 
compromise  between  calculation  time  and  the  requirement  to  convey  the  physics  embodied  in  the 
simulations. 

Evaluation  of  the  equations  requires  the  use  of  the  quad  precision  complex  data  type.  The 
language  of  choice  is  FORTRAN  (version  3.0)  running  on  a  Sun  Microsystems  SPARC  Station 
model  lOSX-512,  Operating  System  Solaris  2.3.  The  calculation  time  varies  from  45  minutes  to 
6  hours. 

The  FORTRAN  algorithms  developed  to  produce  the  simulations  contained  here  are 
published  in  a  separate  document,  which  is  available  via  e-mail  at 

peloquin@cascade.nl.nuwc.navy.mil. 
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RESULTS 

SINGLE-LAYER  CYLINDER  SIMULATIONS— RADIAL  PRESSURE 
EXCITATIONS 

Types  of  Wave  Propagation 

In  review,  an  infinite  isotropic  elastic  media  will  support  two  distinct  types  of  waves:  the 
dilatational  wave  and  the  transverse  (or  shear)  wave.  The  phase  velocity  of  the  dilatational  wave 
is  given  by  equation  (320)  as 


= 


fp  ri+vj(i-2vj- 


(320) 


For  a  dilatational  wave,  the  particle  displacement  and  the  wave  propagation  vector,  which  defines 
the  direction  of  energy  propagation,  are  coincident.  For  the  wavenumber  and  frequency  range 
shown  in  the  simulations,  the  manifestation  of  this  wave  in  the  long  thin  cylinder  produces  a 
modified  phase  velocity,  known  as  the  bar  velocity  or  extensional  wave  phase  velocity.  It  is 
given  in  equation  (321)  as 


c 


e 


(321) 


The  other  wave,  the  transverse  (or  shear)  wave,  is  characterized  by  a  particle  displacement 
that  is  perpendicular  to  the  wave  propagation  vector.  An  expression  for  the  shear  modulus  in 
terms  of  the  Young’s  modulus  and  Poisson  ratio  of  a  material  is  given  by  equation  (322)  as 


= 


E 


C 


2(1 +vj- 


(322) 


From  equation  (11),  the  shear  wave  phase  velocity  is  calculated.  The  phase  velocities  of  the 
dilatational  and  extensional  waves  of  the  cylinder  and  the  dilatational  wave  of  the  fluids  are  listed 
in  table  1  and  are  based  on  the  data  provided  in  table  2. 
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Table  1.  Single-Layer  Cylinder/Fluid  Phase  Velocities 


Wave  Type 

Phase  Velocity 
(m/sec) 

Dilatational  Wave  C/ 

447.51 

Transverse  Wave  c, 

182.68 

Extensional  Wave  Cg 

305.71 

Dilatational  Wave 

1500.0 

The  system  under  consideration  in  this  report  is  infinite  in  longitudinal  coordinate  x;  however, 
it  is  bounded  in  r  and  0.  Furthermore,  the  cylinder  is  in  contact  with  fluids.  This  composite 
system  will  produce  branches  of  wave  propagation  that  are  related  to  the  dilatational  and 
transverse  waves  of  the  media  but  which  do  not  necessarily  propagate  at  the  same  velocity.  The 
velocities  listed  in  table  1  will  serve  as  a  guide  for  analyzing  the  responses  that  follow. 


Mechanical  damping  in  the  cylinder  is  incorporated  by  the  use  of  the  structural  loss  factor  . 
This  structural  loss  factor  represents  the  fraction  of  the  Young’s  modulus  that  causes  the  strain 
response  to  be  out  of  phase  with  the  applied  stress: 


,  (323) 

where  j  equals  1  or  2,  corresponding  to  the  first  and  second  layers  of  the  cylinder, 
respectively.  It  is  therefore  possible,  with  the  formulation  presented  here,  to  have  different 
structural  loss  factors  for  each  cylinder.  This  is,  in  fact,  how  the  two-layer  simulations  were 
performed. 

Cylinder  and  Fluid  Material  Properties 

The  material  properties  for  the  single-layer  cylinder/fluid  simulations  are  listed  in  tables  2 
through  5.  The  values  listed  approximate  water  as  both  an  inner  and  an  outer  fluid.  The 
cylinder  properties  are  those  of  urethane  or  rubber,  much  like  an  ordinary  garden  or  fire  hose. 
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Table  2.  Single-Layer-Cylinder  Properties — ^3.00-in.  Diameter 


Property 

Definition 

E.  =  1.0x10^^ 

Young’s  Modulus 

Cl  =  0.3 

Structural  Loss  Factor 

Pj  =  1070^ 

Density 

II 

o 

Poisson’s  Ratio 

a  =  1.200  in. 

Inner  Radius 

b  =  1.500  in. 

Outer  Radius 

h-^  =  0.300  in. 

Thickness 

Table  3.  Single-Layer-Cylinder  Properties — 0.670-in.  Diameter 


Property 

Definition 

oo 

O 

X 

p 

II 

Young’s  Modulus 

o 

II 

Structural  Loss  Factor 

p,  =  107oiS 

Density 

o 

II 

Poisson’s  Ratio 

a  =  0.235  in. 

Inner  Radius 

b  =  0.335  in. 

Outer  Radius 

/Zj  =  0.100  in. 

Thickness 
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Table  4.  Single-Layer-Cylinder  Outer  Fluid  Properties  (Water) 


Property 

Definition 

=  1000.0^ 

Density 

c,  =  1500.0  ™ 
sec 

Velocity  of  Sound 

Table  5.  Single-Layer-Cylinder  Inner  Fluid  Properties  (Water) 


Property 

Definition 

p.  =  1000.0^ 

Density 

C.  =  1500.0 
'  sec 

Velocity  of  Sound 

Circumferential  Order  Shapes 

For  the  simulations  that  follow,  most  of  the  significant  branches  of  wave  propagation  will 
produce  cylinder  deformation  shapes  as  depicted  in  figures  7  through  10.  Figure  7  depicts  the 
radial  M-displacement  that  occurs  in  the  cylinder  as  energy  is  propagating  in  the  first  branch  of 
wave  propagation,  corresponding  to  circumferential  order  number  n  =  0.  This  first  branch  of 
wave  propagation  is  commonly  referred  to  as  a  breathing  wave.  Figure  8  depicts  the 
longitudinal  w-displacement  that  occurs  in  the  cylinder  as  energy  is  propagating  in  the  second 
branch  of  wave  propagation,  corresponding  to  circumferential  order  number  «  =  0.  This  second 
branch  of  wave  propagation  is  commonly  referred  to  as  an  extensional  wave.  The  displacement 
shapes  depicted  in  figures  9  (n  =  1)  and  10  (n  =  2)  cause  the  cylinder  to  undergo  a  state  of 
bending. 

The  cylinder  deformation  in  figures  7  through  10  has  been  magnified  for  the  purpose  of 
visualizing  the  deformation.  It  should  be  remembered  that  this  is  a  small  strain  analysis,  suitable 
for  analyzing  problems  undergoing  loading  on  an  acoustic  scale. 
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Figure  7.  Circumferential  Order  Number  n  =  0,  First  Branch 


Figure  8.  Circumferential  Order  Number  «  =  0,  Second  Branch 
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Radial  Stress  Variation  in  the  Cylinder  Wall 

This  first  set  of  simulations  is  performed  for  the  data  listed  in  tables  2,  4,  and  5.  The  radial 
stress  transfer  surface  is  calculated  at  three  different  radial  positions  in  figures  11,  12,  and 
13.  These  positions  correspond  to  the  outer,  middle,  and  inner  surface  of  the  cylinder  wall.  In 
reality,  the  inner  and  outer  radial  positions  are  0.01  inch  from  the  actual  free  surface  of  the 
cylinder. 

The  curves  plotted  in  figures  14  and  15  show  the  distribution  and  variation  of  radial  stress  as  a 
function  of  radial  position  r  within  the  cylinder  wall.  The  principal  assumptions  made  in  shell 
theory,^^  which  ignore  the  variations  of  quantities  across  the  cylinder  wall  thickness,  are  not  made 
in  elasticity  theory.  The  model  developed  here  does  not  make  any  assumptions  that  discard  parts 
of  the  physical  behavior  of  the  system.  As  can  be  clearly  observed  in  figures  14  and  15,  the 
radial  stress  varies  across  the  cylinder  wall  thickness. 

In  figure  15,  the  radial  stress  increases  with  decreasing  diameter,  which  amplifies  the  stress 
level  beneath  the  outer  diameter  of  the  cylinder.  This  behavior  is  consistent  with  the  principle  of 
conservation  of  energy  owing  to  the  fact  that  the  circumferential  area  at  a  radial  location  within 
the  cylinder  thickness  is  actually  less  than  the  circumferential  area  at  the  cylinder  outer  diameter. 

The  radial  stress  boundary  conditions  imposed  on  the  cylinder  wall  require  the  stress  to  vary 
from  a  value  equal  to  the  inner  fluid  pressure  field  at  r  =  a  to  a  value  equal  to  the  outer  fluid 
pressure  field  at  r  =  b.  This  variation  can  be  observed  in  the  behavior  of  the  curves  in  figure  14. 

The  other  significant  feature  in  the  simulations  is  the  resonance  at  approximately  2800  Hz, 
which  is  clearly  apparent  in  figure  15.  The  cylinder  is  experiencing  a  predominantly  radial 
vibration  that  is  heavily  influenced  by  the  presence  of  the  fluids  in  this  case.  This  resonance  will 
be  discussed  in  some  detail  in  a  later  section  of  this  report. 
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Figure  11.  Radial  Stress  Transfer  Surface  With  Magnitude  =  10Log(T;./r,)/P(,)^ 
at  rj  =  1.49  in.  for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in. 
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Figure  12.  Radial  Stress  Transfer  Surface  With  Magnitude  =  10Log(t^/r,)/Po)^ 
at  Tj  =  1.35  in.  for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in. 


Figure  13.  Radial  Stress  Transfer  Surface  With  Magnitude  =  10Log(T^/rjI/Po)^ 
at  ri=  1.21  in.  for  Elasticity  Model,  Where  w  =  0  and  Diameter  =  3.00  in. 
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Figure  14.  Comparison  of  Radial  Stress  Transfer  Surfaces  With 
Magnitude  =  10Log(vfri)/P^)2  at  Various  r  for  Elasticity  Model, 
Where  n  =  0,  Diameter  =  3.00  in.,  and/=  200  Hz 
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Figure  15.  Comparison  of  Radial  Stress  Transfer  Surfaces  With 
Magnitude  =  10Log(v(ri)/F„)2  for  Elasticity  Model, 
Where  n  =  0,  Diameter  =  3.00  in.,  and  fc  =  0  rad/m 
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Model  Comparisons — 3.00-Inch-Diameter  Cylinder 

In  this  section,  we  will  compare  the  elasticity  model  developed  in  this  report  to  the  membrane 
shell  and  bending  shell  models  presented  in  reference  13.  As  noted  previously,  the  assumptions 
made  for  shell  theory  have  limitations  on  the  mean  diameter-to-thickness  ratio  (which  should  be 
above  10  for  accurate  model  behavior)  as  well  as  on  the  wavelength-to-diameter  and  wavelength- 
to-thickness  ratios.  The  elasticity  model  represents  the  true  behavior  of  the  cylinder  because  it 
does  not  introduce  simplifying  assumptions  for  the  elastic  behavior  of  the  material.  Although  it 
is  acknowledged  that  the  shell  models  will  not  always  be  in  full  agreement  with  the  elasticity 
model,  a  further  understanding  of  the  fluid/structure  interaction  is  obtained  by  investigating  the 
agreement  that  does  occur  between  the  models  for  the  3.00-inch-outer-diameter  cylinder  (in  this 
section)  and  the  0.670-inch-outer-diameter  cylinder  (in  the  following  section). 

Each  simulation  is  displayed  as  a  color  image,  where  color  indicates  magnitude.  For  clarity, 
and  in  an  effort  to  reveal  the  true  character  of  the  surface,  the  color  image  is  also  displayed  as  a 
wire  frame  surface  directly  below  the  corresponding  color  image. 

Figures  16  through  21  are  the  simulations  {n  =  0)  for  the  membrane  shell,  bending  shell,  and 
elasticity  models,  respectively.  In  figure  22,  which  is  a  comparison  of  the  three  models  at  ^  =  0, 
both  shell  models  predict  the  resonance  frequency  of  the  cylinder  to  be  higher  than  does  the 
elasticity  model.  The  shell  models  indicate  the  resonance  frequency  to  be  at  approximately  11 
kHz,  while  the  elasticity  model  is  indicating  the  resonance  at  2870  Hz.  The  combination  of 
urethanelike  properties  (cylinder)  and  water  (fluids)  causes  the  shell  models  to  incorrectly  predict 
the  first  radial  resonance  of  the  cylinder.  If  air  had  been  chosen  for  the  inner  and  outer  fluids,  the 
resonance  frequencies  of  the  models  would  have  been  in  close  agreement. 

The  breathing  wave,  evident  in  the  three  model  types,  is  dispersive  in  all.  Figures  23  and  24 
compare  the  breathing  wave  resonance  peaks  at  200  and  600  Hz,  respectively.  The 
corresponding  phase  velocity  for  these  two  frequencies  is  94  and  74  m/sec  for  the  elasticity 
model,  and  126  and  80  m/sec  for  the  shell  model.  At  high  frequency,  the  character  of  the 
breathing  wave  in  the  membrane  model  is  very  different  from  both  the  bending  shell  and  elasticity 
models.  Above  1000  Hz,  the  breathing  wave  in  the  membrane  model  slows  dramatically  and 
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falls  in  amplitude  relative  to  the  other  two  models.  In  figure  23,  at  80  rad/m,  the  membrane 
model  begins  to  underpredict  the  transfer  function  roll-off,  ultimately  underpredicting  the  roll-off 
relative  to  the  elasticity  model  prediction  by  approximately  40  dB  at  300  rad/m. 

The  extensional  wave  is  evident  in  the  elasticity  model  of  figure  24.  It  appears  as  low-level 
peaks  on  either  side  of  A:  =  0  at  a  phase  velocity  of  291  m/sec  and  exhibits  no  dispersion.  In 
figure  25,  the  extensional  wave  is  more  clearly  visible  at  100  rad/m;  it  is  also  visible  in  figure 
20.  Comparison  of  the  three  models  in  figure  25  indicates  significant  divergence  of  the  shell 
solutions  from  the  elasticity  solutions  above  20  rad/m. 

The  simulations  for  circumferential  order  number  «  =  1  are  shown  in  figures  26  through 
31.  The  first  branch  of  wave  propagation  places  the  cylinder  in  a  state  of  bending,  as  shown  in 
the  wire  frame  representation  of  figure  9.  Two  additional  branches  of  wave  propagation  are 
visible  in  the  simulations.  The  second  branch  has  a  cutoff  frequency  of  approximately  1000  Hz 
and  causes  the  longitudinal  displacement  to  undergo  one  wavelength  of  variation  around  the 
circumference  of  the  cylinder.  The  third  branch  has  a  cutoff  frequency  of  approximately  2000 
Hz,  causing  the  circumferential  displacement  to  undergo  one  wavelength  of  variation  around  the 
circumference  of  the  cylinder.  Both  of  these  branches  couple  into  the  radial  displacement, 
which  develops  a  pressure  in  the  inner  fluid.  However,  these  pressure  field  components  are 
weak  because  of  the  magnitude  of  the  damping  =  0.3).  The  reader  is  referred  to  reference 
13  for  a  complete  set  of  figures  depicting  the  deformation  shapes  of  the  cylinder  for  these 
branches  of  wave  propagation.  Figures  32  through  34  compare  the  three  models  dXk-0  and  at 
various  cuts  in  wavenumber.  The  resonant  peak  of  the  first  branch  of  wave  propagation 
compares  favorably  for  all  three  models  up  to  600  Hz,  as  seen  in  figure  34.  The  membrane 
model’s  roll-off  in  wavenumber  departs  from  the  other  two  models  by  approximately  60  rad/m, 
and  at  4500  Hz  (figure  35)  the  significant  divergence  of  the  shell  models  from  the  elasticity  model 
renders  the  shell  models  to  be  unacceptable  predictors  of  the  cylinder’s  behavior. 

Figures  36  through  41  are  the  simulations  for  circumferential  order  number  n  =  2.  We  see 
from  figures  36  and  42  that  the  membrane  model  will  not  support  a  cutoff  frequency  for  the  first 
branch  of  wave  propagation.  This  first  branch  produces  a  deformation  in  the  cylinder  as  shown 
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in  figure  10.  Two  additional  branches  are  visible,  both  with  cutoff  frequencies.  Each  produces 
the  same  kind  of  deformation  in  the  shell  as  described  for  n  =  1,  only  now  there  are  two 
wavelengths  of  variation  around  the  circumference.  In  figures  43  through  45,  we  see  that  the 
bending  shell  model  is  useful  up  to  600  Hz,  diverging  from  the  elasticity  model  only  at  higher 
frequencies. 

To  summarize,  for  the  3.00-inch-diameter  cylinder  considered  here,  the  shell  models  are  in 
good  agreement  with  the  elasticity  model  for  the  region  extending  to  1000  Hz  and  ±75  rad/ 
m.  Outside  of  this  region,  the  difference  increases  beyond  5  dB,  and  the  elasticity  solution 
should  be  used.  For  the  n  =  0  case,  the  fluid  loading  influences  the  first  radial  natural  frequency 
of  the  shell  models  so  dramatically  that  the  error  relative  to  the  elasticity  model  is  greater  than  a 
factor  of  three,  as  evidenced  by  figure  22.  For  the  n  =  2  case,  the  membrane  model  should  not  be 
used,  owing  to  the  misrepresentation  of  the  behavior  of  the  first  branch  (i.e.,  lack  of  cutoff 
frequency  and  appropriate  branch  structure).  Although  not  shown  here,  similar  behavior  will 
occur  with  the  membrane  shell  model  for  n  >  2 . 
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Figure  16.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/ri)/P^)2  at  -  0.75  in. 
for  Membrane  Shell  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  17.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)^  at  =  0.75  in. 
for  Membrane  Shell  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  18.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/r,)/P^)2  at  r,  =  0.75  in. 
for  Bending  Shell  Model,  Where  «  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  19.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P„)2  at  rj  =  0.75  in. 
for  Bending  Shell  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  20.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P^)^  at  r,  =  0.75  in. 
for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  21.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(ri)/P^)^  at  r|  =  0.75  in. 
for  Elasticity  Model,  Where  «  =  0  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  22.  Model  Comparison  for  Figures  16, 18,  and  20  at  Extended  Frequency 
Range  When  Magnitude  =  10Log(P/ri)/P<,)2  at  rj  =  0.75  in.  With  n  =  0, 
Diameter  =  3.00  in.,  and  ^  =  0  rad/m 
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Figure  23.  Model  Comparison  for  Figures  16, 18,  and  20  When 
Magnitude  =  10Log(P/ri)/P^)^  at  rj  =  0.75  in.  With  n  =  0, 
Diameter  =  3.00  in.,  and  /=  200  Hz 
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Figure  24.  Model  Comparison  for  Figures  16, 18,  and  20  When 
Magnitude  =  10Log(P/rij/F«)2  at  =  0.75  in.  With  n  =  0, 
Diameter  =  3.00  in.,  and  /=  600  Hz 
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Figure  25.  Model  Comparison  for  Figures  16, 18,  and  20  When 
Magnitude  =  lOLogCF/rJ/F^)^  at  =  0.75  in.  With  #i  =  0, 
Diameter  =  3.00  in.,  and  /  =  4500  Hz 


IS 


WAVENUMBER  (rad/ 


TR  1 1 ,067 


200.00 


-200.00 


MAGNITUDE 


1000.00  2000.00  3000.00 

FREQUENCY  (Hz) 


4000.00  5000.00 


Figure  26.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r|)/P^)^  at  r^  -  0.75  in. 
for  Membrane  Shell  Model,  Where  n  =  1  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  27.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P(,)^  at  r^  =  0.75  in. 
for  Membrane  Shell  Model,  Where  n-  \  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  28.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P^)^  at  rj  =  0.75  in. 
for  Bending  Shell  Model,  Where  «  =  1  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  29.  Pressure  Transfer  Surface  With  Magnitude  =  lOLogCF/rjj/P^)^  at  rj  =  0.75  in. 
for  Bending  Shell  Model,  Where  n  =  1  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  30.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/ri)/P„)2  at  rj  =  0.75  in. 
for  Elasticity  Model,  Where  n  =  1  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  31.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P,<ri)/Po)^  at  '‘i  =  0.75  in. 
for  Elasticity  Model,  Where  n  =  1  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  32.  Model  Comparison  for  Figures  26, 28,  and  30  When 
Magnitude  =  10Log(P,-(ri)/Po)^  at  =  0.75  in.  With  n  =  1, 
Diameter  =  3.00  in.,  and  A:  =  0  rad/m 
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Figure  34.  Model  Comparison  for  Figures  26, 28,  and  30  When 
Magnitude  =  lOLogCF/rj)//*^)^  at  =  0.75  in.  With  n  =  1, 
Diameter  =  3.00  in.,  and/  =  600  Hz 
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Figure  36.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/(ri)/Po)^  at  rj  =  0.75  in. 
for  Membrane  Shell  Model,  Where  n  =  2  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  37.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P,-(ri)/P„)^  at  rj  =  0.75  in. 
for  Membrane  Shell  Model,  Where  n  =  2  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  38.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)2  at  r,  =  0.75  in. 
for  Bending  Shell  Model,  Where  «  =  2  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  39.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/Po)2  at  ri  =  0.75  in. 
for  Bending  Shell  Model,  Where  «  =  2  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  40.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/r at  rj  -  0.75  in. 
for  Elasticity  Model,  Where  n  =  2  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  41.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/rD/P<,)^  at  -  0.75  in. 
for  Elasticity  Model,  Where  n  =  2  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  42.  Model  Comparison  for  Figures  36, 38,  and  40  When 
Magnitude  =  10Log(P,(r|>/Po)^  at  ri  =  0.75  in.  With  n  =  2, 
Diameter  =  3.00  in.,  and  A:  =  0  rad/m 
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Figure  43.  Model  Comparison  for  Figures  36, 38,  and  40  When 
Magnitude  =  10Log(Fj('rj)/P^)^  at  =  0.75  in.  With  «  =  2, 
Diameter  =  3.00  in.,  and  /  =  200  Hz 
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Model  Comparisons — 0.670-Inch-Diameter  Cylinder 

Let  us  now  consider  the  case  of  a  cylinder  with  the  properties  listed  in  table  3,  surrounded  by 
fluids  with  the  properties  listed  in  tables  4  and  5.  This  cylinder,  similar  to  the  one  previously 
considered,  has  properties  that  are  representative  of  a  garden  or  Are  hose,  and  the  fluids  have 
properties  consistent  with  those  of  water.  Figures  46  through  5 1  are  the  n  =  0  simulations  for  the 
three  models.  The  output  quantity  is  the  pressure  in  the  inner  fluid;  this  field  quantity  is 
calculated  at  rj  =  0.1  inch.  Figures  52  through  55  are  the  comparisons  of  the  three 
models.  The  major  branch  of  wave  propagation  is  the  breathing  wave,  indicated  on  figures  47, 
49,  and  51.  The  extensional  wave  is  visible  in  the  elasticity  model  (figures  51  and  55).  All 
three  models  are  in  close  agreement  up  to  ±150  rad/m.  Above  this  point,  the  membrane  shell 
model  diverges;  however,  the  bending  shell  model  continues  to  remain  within  4  dB  of  the 
elasticity  model  up  to  4500  Hz  and  ±300  rad/m. 

The  large  discrepancy  between  the  shell  models  and  the  elasticity  model  for  the  first  radial 
natural  frequency  of  the  cylinder  is  apparent  in  figure  52.  The  fluid-loaded  shell  models  produce 
resonances  at  53,136  Hz,  while  the  fluid-loaded  elasticity  model  displays  this  resonance  at  10,793 
Hz.  As  was  seen  in  the  3.00-inch-diameter  simulations  of  the  previous  section,  the  shell  models 
are  not  adequate  for  predicting  the  first  radial  natural  frequency  when  the  present  combination  of 
water  and  rubberlike  cylinder  properties  is  chosen.  The  elasticity  model  provides  the  best 
solution  in  this  case. 

Figures  56  through  61  are  the  simulations  for  circumferential  order  number  n  =  1.  Figure  62 
is  a  comparison  of  the  models  across  frequency,  and  figures  63  through  65  are  comparisons  of  the 
models  across  wavenumber  for  particular  frequencies.  The  first  and  dominant  branch  of  wave 
propagation  in  the  simulations  is  one  in  which  the  cylinder  is  placed  in  a  state  of  bending,  as 
depicted  by  the  wire  frame  representation  of  figure  9.  The  shell  models  are  in  close  agreement 
with  the  elasticity  model  up  to  ±150  rad/m  and  5000  Hz,  as  seen  in  figures  62  through 
65.  Above  these  points  in  the  wavenumber-frequency  plane,  the  membrane  shell  model 
diverges.  The  bending  shell  model,  however,  remains  within  4  dB  of  the  elasticity  solution  over 
the  entire  range  of  wavenumbers  and  frequencies  considered. 
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The  next  branch  of  wave  propagation,  faintly  visible  in  the  4500-Hz  region,  has  a  cutoff 
frequency  associated  with  it.  As  described  earlier  for  the  3.00-inch  cylinder,  this  second  branch 
is  causing  the  longitudinal  displacement  to  undergo  one  wavelength  of  variation  with  respect  to 
the  circumference  of  the  cylinder.  The  predominantly  longitudinal  displacement  does  couple 
into  the  radial  displacement,  thereby  causing  an  effect  on  the  pressure  field  in  the  inner 
fluid.  The  value  of  the  damping  selected  for  the  cylinder  has  diminished  the  sharpness  of  this 
branch. 

Figures  66  through  71  contain  the  simulations  for  circumferential  order  number  n  =  2.  As  in 
the  case  for  the  cylinder  considered  in  the  previous  section,  the  membrane  shell  model  will  not 
support  a  cutoff  frequency  for  the  first  branch  of  wave  propagation  (figures  66  and  72).  This  is 
due  to  the  assumption  that  neglects  bending  stiffness  terms,  which  are  inherent  in  a  membrane 
formulation.  Figures  73  through  75  are  comparisons  across  wavenumber,  and  once  again  we  see 
the  divergence  of  the  membrane  shell  model  and  the  close  agreement  between  the  bending  shell 
and  elasticity  solutions. 

To  summarize,  the  membrane  shell  model  generally  falls  within  4  dB  of  the  elasticity  model 
for  the  n  =  0  circumferential  order  number  in  the  region  +150  rad/m  and  0  to  5000  Hz.  If  it  is 
necessary  to  obtain  a  simulation  with  an  error  of  less  than  4  dB,  the  elasticity  model  should  be 
used.  Furthermore,  the  elasticity  model  (not  the  shell  models)  should  also  be  used  for  the 
determination  of  the  first  radial  natural  frequency  of  this  cylinder/fluid  combination.  For  the 
higher  circumferential  order  numbers,  it  should  be  remembered  that  the  membrane  model  will  not 
support  a  cutoff  frequency  for  the  first  branch  of  wave  propagation.  The  bending  shell  model 
will  be  within  4  dB  of  the  elasticity  solution.  If  an  error  of  less  than  4  dB  is  required,  the 
elasticity  model  should  be  used. 
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Figure  46.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/F^)2  at  rj  =  0.1  in. 
for  Membrane  Shell  Model,  Where  n  =  0  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  47.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(rD/P<,)^  at  U  =  0-1  in. 
for  Membrane  Shell  Model,  Where  n  =  0  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  48.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P,(ri)/Po)^  G  =  0.1  in. 
for  Bending  Shell  Model,  Where  n  =  0  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  49.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P(,)^  G  =  0.1  in. 
for  Bending  Shell  Model,  Where  n  =  0  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  50.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)^  at  r,  =  0.1  in. 
for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  51.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj<r,)/P^)^  at  r,  =  0.1  in. 
for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  52.  Model  Comparison  for  Figures  46, 48,  and  50  at  Extended  Frequency 
Range  When  Magnitude  =  10Log(Pj(r,)/F„)^  at  =  0.1  in.  With  n  =  0, 
Diameter  =  0.670  in.,  and  A:  =  0  rad/m 
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Figure  54.  Model  Comparison  for  Figures  46, 48,  and  50  When 
Magnitude  =  10Log(P/ri)/F^)2  at  =  0.1  in.  With  u  =  0, 
Diameter  =  0.670  in.,  and/=  600  Hz 
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Figure  56.  Pressure  Transfer  Surface  With  Magnitude  =  lOLogCP/rD/P,,)^  at  =  0.1  in. 
for  Membrane  Shell  Model,  Where  n  =  1  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  57.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P,(ri)/P<,)^  at  U=  0-1  in. 
for  Membrane  Shell  Model,  Where  n  =  1  and  Diameter  0.670  in.  (Wire  Frame  Surface) 


100 


WAVENUMBER  (rad/ 


TR  11,067 


Figure  58.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/r,)/P^)^  at  =  0.1  in. 
for  Bending  Shell  Model,  Where  «  =  1  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  59.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)^  at  =  0.1  in. 
for  Bending  Shell  Model,  Where  n  =  1  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  60.  Pressure  Transfer  Surface  With  Magnitude  =  10  Log(P/ri)/F^)^  at  rj  =  0.1  in. 
for  Elasticity  Model,  Where  n  =  1  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  61.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Fj(ri)/PJ^  at  =  0.1  in. 
for  Elasticity  Model,  Where  n-  \  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  62.  Model  Comparison  for  Figures  56, 58,  and  60  When 
Magnitude  =  10Log(P,(ri)/P^)2  at  =  0.1  in.  With  n  =  1, 
Diameter  =  0.670  in.,  and  k  =  0  rad/m 
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Figure  63.  Model  Comparison  for  Figures  56, 58,  and  60  When 
Magnitude  =  10Log(Pjfri)/P^)^  at  rj  =  0.1  in.  With  n  =  1, 
Diameter  =  0.670  in.,  and  /=  200  Hz 
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Figure  65.  Model  Comparison  for  Figures  56, 58,  and  60  When 
Magnitude  =  10Log(P/ri)/F^)^  at  =  0.1  in.  With  «  =  1, 
Diameter  =  0.670  in.,  and/=  4500  Hz 


106 


WAVENUMBER  (rad/m 


TR  11,067 


MAGNITUDE 


600.00 

400.00 

200.00 

0.00 

-200.00 

-400.00 

-600.00 


-12 

-15 


1000.00  2000.00  3000.00  4000.00  5000.00 

FREQUENCY  (Hz) 


Figure  66.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P^)^  at  T]  =  0.1  in. 
for  Membrane  Shell  Model,  Where  n  =  2  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  67.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(r|)/P^)^  at  r^=  0.1  in. 
for  Membrane  Shell  Model,  Where  n  =  2  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  68.  Pressure  Transfer  Surface  With  Magnitude  =  lOLogCP/ri/P^)^  at  =  0.1  in. 
for  Bending  Shell  Model,  Where  n  =  2  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  69.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P^)^  at  =  0.1  in. 
for  Bending  Shell  Model,  Where  n  =  2  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  70.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)^  at  r,  =  0.1  iu. 
for  Elasticity  Model,  Where  w  =  2  and  Diameter  =  0.670  in.  (Color  Image) 


Figure  71.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(r|)/P^)^  at  =  0.1  in. 
for  Elasticity  Model,  Where  n-2  and  Diameter  =  0.670  in.  (Wire  Frame  Surface) 
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Figure  72.  Model  Comparison  for  Figures  66,  68,  and  70  When 
Magnitude  =  10Log(Pj(r^)/P^)^  at  =  0.1  in.  With  n  =  2, 
Diameter  =  0.670  in.,  and  fc  =  0  rad/m 
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Figure  74.  Model  Comparison  for  Figures  66, 68,  and  70  When 
Magnitude  =  10Log(F,<ri)/P^)^  at  =  0.1  in.  With  n  =  2, 
Diameter  =  0.670  in.,  and/ =  600  Hz 
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Inner  Fluid  Pressure  Field  Evanescence 

From  this  point  on,  unless  otherwise  noted,  the  elasticity  model  will  be  used  for  the 
simulations  conducted  in  the  exploration  of  the  cylinder/fluid  dynamic  behavior. 

Let  us  consider  the  effect  of  the  propagating  region  and  the  nonpropagating  region  on  the 
magnitude  of  the  pressure  field  that  exists  in  the  inner  fluid  as  a  function  of  radial  position.  For 
wavenumbers  less  than  co/Cp  the  pressure  field  propagates  as  shown  by  equation  (67).  For 
wavenumbers  greater  than  co/Cj-,  the  pressure  field  becomes  increasingly  evanescent,  with 
increasing  wavenumber  and  decreasing  radial  position  r,  as  shown  by  equation  (71). 

The  simulation  for  the  inner  fluid  pressure  field  at  rj  =  1.1  inches  is  shown  in  figures  76  and 
77.  The  simulation  of  this  pressure  field  at  =  0.01  inch  is  shown  in  figures  78  and  79.  In 
figures  80,  81,  and  82,  comparison  of  the  pressure  field  is  made  for  these  two  radial  positions  at 
three  different  frequencies.  Attenuation  of  the  pressure  field,  which  is  apparent  in  all  the  figures, 
is  on  the  order  of  12  dB  at  1(K)  rad/m. 

Comparison  of  the  images  in  figures  76  and  78  reveals  a  decrease  in  breathing  wave  resonant 
peak  level  with  increasing  wavenumber.  The  breathing  wave  pressure  is  locally  generated  at  the 
cylinder/fluid  boundary.  Because  the  breathing  wave  phase  velocity  is  an  order  of  magnitude 
slower  than  the  inner  fluid  propagation  velocity,  the  magnitude  of  the  pressure  will  suffer 
increasing  attenuation  away  from  the  surface  of  generation  as  the  breathing  wave  diverges  in 
wavenumber  (compare  figures  80  and  81). 

Similar  behavior  is  observed  with  the  extensional  wave.  However,  since  the  phase  velocity 
of  the  extensional  wave  is  greater  than  the  breathing  wave,  similar  attenuation  does  not  occur 
until  a  higher  frequency  is  reached.  In  figure  82  (4500  Hz),  the  extensional  wave  has  undergone 
approximately  the  same  attenuation  (12  dB)  as  the  breathing  wave  has  undergone  in  figure  81  (12 
dB)  at  900  Hz. 
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Figure  76.  Pressure  Transfer  Surface  With  Maguitude  =  10Log(P/r,)/P^)^  at  r,  =  1.1  in. 
Where  «  =  0  and  Diameter  =  3.00  in.  for  Field  Decay  Comparison  (Color  Image) 


Figure  77.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/Po)^  at  U  =  1-1  in. 
Where  n  =  0  and  Diameter  =  3.00  in.  for  Field  Decay  Comparison  (Wire  Frame  Surface) 
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Figure  78.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/rj)/P at  =  0.01  in. 
Where  n  =  0  and  Diameter  =  3.00  in.  for  Field  Decay  Comparison  (Color  Image) 


Figure  79.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F,(ri)/Po)^  at  =  0.01  in 
Where  «  =  0  and  Diameter  =  3.00  in.  for  Field  Decay  Comparison  (Wire  Frame  Surface) 
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Figure  81.  Field  Decay  Comparison  for  Figures  76  and  78  When 
Magnitude  =  10Log(P,<r|)/P<,)^  at  Various  rj  With  n  =  0, 
Diameter  =  3.00  in.,  and  /  =  900  Hz 


Figure  82.  Field  Decay  Comparison  for  Figures  76  and  78  When 
Magnitude  =  10Log(P,fri)/P^)^  at  Various  rj  With  n  =  0, 
Diameter  =  3.00  in.,  and  /  =  4500  Hz 
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Inner  Fluid  Density  Variation 

In  this  section,  we  shall  explore  the  effects  of  a  change  in  inner  fluid  density  on  the  attenuation 
of  the  inner  fluid  pressure  field.  The  plane  wave  velocity  or  dilatational  velocity  of  the  inner 
fluid  will  remain  fixed  at  1500  m/sec.  Simulations  performed  with  an  inner  fluid  density  of  500, 
1000,  and  2000  kg/m^  are  displayed  in  figures  83,  84  and  85. 

In  figure  86,  comparisons  are  made  at  ^  =  0  rad/m.  The  decrease  in  inner  fluid  density  has 
lowered  the  first  radial  resonance  of  the  cylinder/fluid  combination.  Figures  87  through  89  are 
comparisons  of  the  effect  of  the  change  in  density  across  wavenumber  at  three  different 
frequencies.  In  figure  87,  we  see  that  a  decrease  of  inner  fluid  density  has  increased  the  phase 
velocity  of  the  breathing  wave.  Breathing  wave  phase  velocity  from  figure  87  at  300  Hz  is 
provided  in  table  6. 

Table  6.  Effect  of  Inner  Fluid  Density  on  Breathing  Wave  Phase  Velocity  From  Figure  87 


Inner  Fluid  Density 
(kg/m^) 

Breathing  Wave 

Phase  Velocity 
(m/sec) 

500 

130 

1000 

91 

2000 

66 
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Figure  83.  Pressure  Transfer  Surface  With  Magnitude  =  lOLogCP/rD/P^)^  at  r,  =  0.01  in. 
for  Pi  =  500  kg/m-^,  n  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 
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Figure  84.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)2  at  r,  =  0.01  in. 
for  Pi  =  1000  kg/m^,  n  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 
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Figure  85.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r|)/P^)^  at  rj  =  0.01  in, 
for  pj  =  2000  kg/m^,  «  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 
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Figure  86.  Effect  of  Inner  Fluid  Density  When  Magnitude  =  lQLog{P 
at  1*1  =  0.01  in.  for  n  =  0,  Diameter  =  3.00  in.,  and  =  0  rad/m 
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Figure  89.  Effect  of  Inner  Fluid  Density  When  Magnitude  =  10Log(P,(ri)/Po)^  at 
Ty  =  0.01  in.  for  n  =  0,  Diameter  =  3.00  in.,  and/  =  4500  Hz 


TR  11,067 


High-Frequency  Response 

Elasticity  and  Bending  Shell  Model  Comparison.  In  previous  sections,  we  have  compared 
the  shell  models  with  the  elasticity  model  up  to  5000  Hz,  accompanied  by  a  limited  comparison 
up  to  15  kHz  at  A:  =  0.  In  this  section,  we  will  compare  the  bending  shell  model  (figures  90  and 
91)  with  the  elasticity  model  (figures  92  and  93)  over  the  range  ±200  rad/m  and  10  to  60,000 
Hz.  These  simulations  have  been  performed  with  the  data  listed  in  tables  2,  4,  and  5. 

The  first  radial  resonance  can  be  clearly  observed  in  the  ^  =  0  cut  displayed  in  figure  94.  As 
was  observed  previously  in  figure  22,  the  bending  shell  model  considerably  overpredicts  the  first 
radial  resonance  of  the  cylinder  because  of  the  fluid  loading.  Both  models  exhibit  deep  nulls  at 
20,704  and  47,214  Hz,  which  are  due  to  the  nature  of  the  inner  fluid  field,  the  propagation 
velocity  of  the  inner  fluid,  and  the  dimensions  of  the  cylinder.  This  behavior  will  be  examined  in 
the  next  section. 

The  remaining  resonances  seen  in  figure  94  (elasticity  model)  occur  as  a  dilatational  wave 
propagates  through  the  thickness  of  the  cylinder  wall.  These  resonances  correspond  to  the  1/2, 
1,  and  3/2  wavelengths  existing  across  the  cylinder  wall  thickness.  This  will  be  explored  in  a 
later  section  for  an  air-loaded  (rather  than  a  water-loaded)  cylinder  with  a  variation  of  wall 
thickness  and  Young’s  modulus. 

The  assumptions  made  in  shell  theory  result  in  an  inadequate  model  for  investigating  the  high- 
frequency  response  of  the  cylinder/fluid  combination  chosen  here.  For  the  heavy  fluid  loading 
present  in  this  combination,  it  is  necessary  to  use  the  elasticity  formulation  to  properly  model  the 
dynamics  of  the  system,  even  at  3000  Hz,  which  is  only  moderately  high  frequency. 
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Figure  90.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(Fj(ri)/P^)^  at 
ri  =  1.1  in.  for  Bending  Shell  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 
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Figure  91.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(P/rD/P(,)^  at 
rj  =  1.1  in.  for  Bending  Shell  Model,  Where  n  -  0  and  Diameter  =  3.00  in. 

(Wire  Frame  Surface) 
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Figure  92.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(P/ri)/F^)^  at 
Vi  =  1.1  in.  for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  93.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(P/ri)/Pj2  at 
=  1.1  in.  for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in. 

(Wire  Frame  Surface) 
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Figure  94.  Model  Comparison  Up  to  60  kHz  for  Figures  90  and  92 
When  Magnitude  =  10Log(P/ri)//*^)^  at  =  1.1  in.  With  n  =  0, 
Diameter  =  3.00  in.,  and  =  0  rad/m 
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Inner  Fluid  Pressure  Field  Nulls.  In  this  section,  we  shall  consider  the  cause  of  the  nulls 
evident  in  figure  94.  For  «  =  0  and 

2 


(the  case  of  propagating  wavenumbers),  the  inner  fluid  field  is  defined  in  terms  of  the  Bessel 
function  as  given  by  equation  (66).  The  nulls  in  frequency  observed  in  figure  94  correspond  to 
the  roots  of  J ^  .  The  argument  of  , 


/• 


1 


(324) 


is  rearranged  in  order  to  solve  for  the  frequency/,,  corresponding  to  the  roots  of  at  yf  =  0 
and  r,  as  follows: 


(325) 


Table  7  lists  the  frequencies  f„  (calculated  with  equation  (325))  that  correspond  to  the  first  two 
roots  of  for  three  different  values  of  radial  position  /•].  A  simulation  of  the  inner  fluid  pressure 


Table  7.  Null  Frequencies  Corresponding  to  the  First  Two  Roots  of 


Root  (n) 

n 

/«at 

Tj  =  1.1  in. 

/«at 

7'i  =  0.55  in. 

/«at 

/■j  =  0.01  in. 

(Hz) 

(Hz) 

(Hz) 

1 

2.4048 

20,545 

(20,704) 

41,092 

(41,056) 

2,259,000 

2 

5.5201 

47,161 

(47,214) 

^  94,326 

5,187,000 

131 


TR  11,067 


field  is  shown  in  figures  95  and  96  evaluated  at  r,  =  0.55  inch.  In  a  similar  manner,  the  pressure 
field  in  the  inner  fluid  at  =  0.01  inch  is  computed  and  displayed  in  figures  97  and 
98.  Comparison  of  the  three  radial  positions  from  figures  92,  95,  and  97  is  made  in  figure  99  at  k 
=  0.  In  figure  99,  the  very  dramatic  effect  of  radial  location  on  the  pressure  field  in  the  inner 
fluid  can  be  observed.  The  nulls  listed  in  table  7  for  =  0.01  inch  occur  at  frequencies  that 
extend  beyond  the  60-kHz  range  of  the  figure  and  are  therefore  unobservable.  The  three  nulls 
observable  in  figure  99  are  listed  in  table  7  in  parentheses  beneath  the  value  calculated  from 
equation  (325).  The  agreement  is  within  1  percent. 

The  null  frequencies/,,  (equation  (325))  will  decrease  if  the  inner  fluid  propagation  velocity  is 
decreased.  Such  will  be  the  case  in  the  next  section,  where  air  is  substituted  for  water  to 
decrease  the  effect  of  the  fluid  on  the  structure  in  order  to  explore  the  second,  third,  and  fourth 
resonances  in  figure  94. 
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Figure  95.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  lOLogCF/rD/F^,)^ 
at  Tj  =  0.55  in.  With  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  96.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(P/rij/Po)^ 
at  Tj  =  0.55  in.  With  n  =  0  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  97.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(Pj(ri)/Po)^ 
at  rj  =  0.01  in.  With  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  98.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(P,(ri)/Po) 
at  r.  =  0.01  in.  With  n  =  0  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  99.  Comparison  of  Inner  Fluid  Pressure  Field  for  Figures  92, 95,  and  97 
When  Magnitude  =  10Log(F,-(ri)/Pp)^  at  Various  With  n  =  0, 
Diameter  =  3.00  in.,  and  k  =  0  rad/m 
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Cylinder  With  Inner  and  Outer  Fluids  of  Air — Wall  Thickness  Variation.  In  this 
section,  we  will  examine  the  three  resonances  above  the  first  radial  resonance  of  the  cylinder,  as 
seen  in  figure  99.  The  wall  thickness  will  be  varied  and  the  effect  of  this  variation  on  the 
resonances  will  be  examined.  The  simulations  in  this  section  and  the  one  immediately  following 
will  be  performed  with  a  change  in  the  inner  and  outer  fluid  properties  from  water  to  air,  as  shown 
in  tables  8  and  9.  As  noted  earlier,  air  was  chosen  in  order  to  reduce  the  fluid’s  effect  on  the 
cylinder.  Although  sharp  nulls  and  peaks  occur  at  lower  frequency  due  to  the  decrease  in  c,- 
(discussed  in  the  last  section),  this  behavior  is  not  strong  enough  to  affect  the  location  of  the 
resonances  under  consideration. 


Table  8.  Outer  Fluid  Properties  for  Wall  Thickness  Variation  (Air) 


Property 

Definition 

P.  =  1-2^ 

Density 

c,  =  343.0  ^ 

^  sec 

Velocity  of  Sound 

Table  9.  Inner  Fluid  Properties  for  Wall  Thickness  Variation  (Air) 


Property 

Definition 

P,-  =  1-2^^ 

Density 

C.  =  343.0  ^ 

‘  sec 

Velocity  of  Sound 

Simulations  for  the  radial  stress  are  shown  in  figures  100  and  101  for  cylinders  of  two 
different  wall  thicknesses.  A  comparison  of  these  two  figures  is  shown  in  figure  102  at 
k  =  0.  The  wall  thickness  was  increased  from  0.3  to  0.6  inch  in  figure  101  so  that  additional 
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resonances  would  be  present  in  the  60-kHz  frequency  range.  From  table  1,  we  see  that  the 
dilatational  phase  velocity  for  the  cylinder  material  is  Cj  =  447.51  m/s.  For  each  resonance 
frequency  observed  in  figures  100  and  101,  the  corresponding  wavelength  L  is  calculated  using 


and  is  provided  in  tables  10  and  1 1 .  The  resonances  correspond  to  dilatational  wave  propagation 
through  the  thickness  of  the  cylinder,  with  standing  waves  existing  in  the  thickness 
direction.  Mode  shapes  of  L/2,  L,  and  3L/2  correspond  to  the  three  resonances  of  figure  101, 
while  only  the  L/2  mode  shape  corresponds  to  the  resonance  in  figure  100.  Doubling  the 
thickness  of  the  cylinder  doubles  the  wavelength  of  the  first  mode  and  halves  the  frequency, 
which  is  observed  with  the  first  resonance,  as  shown  in  tables  10  and  1 1 . 


Table  10.  Thickness  Resonance  Table  for  hy  =  0.3  in.  From  Figure  100 


Thickness 

fti 

Resonance 

L 

hJL 

hy/L  Nominal 

Frequency 

(Hz) 

(m) 

Number  (/) 

1 

30,264 

0.01478 

0.52 

1 

2 

Table  11.  Thickness  Resonance  Table  for  hy  =  0.6  in.  From  Figure  101 


Thickness 
Resonance 
Frequency 
Number  (/) 

fti 

(Hz) 

L 

(m) 

hylL 

hy/L  Nominal 

1 

15,198 

0.02940 

0.52 

1 

2 

2 

30,532 

0.01465 

1.04 

1 

3 

45,221 

0.00989 

1.54 

3 

2 
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Figure  100.  Radial  Stress  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  lOLogCx^^^^frjj/P^)^  at  =  1.21  in.  With  =  0.3  in.,  n  =  0, 
and  Diameter  =  3.00  in.  (Color  Image) 


Figure  101.  Radial  Stress  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  10Log(Xy^^^(ri)/P^)^  at  =  1.21  in.  With  hy  =  0.6  in.,  n  =  0, 
and  Diameter  =  3.00  in.  (Color  Image) 
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Figure  102.  Effect  of  Variation  of  Wall  Thickness  on  the  Cylinder  Response  for 
Figures  100  and  101  When  Magnitude  =  10Log(x^^^^(ri)/F^)^  at  rj  =  1.21  in. 
With  #1  =  0,  Diameter  =  3.00  in.,  and  ^  =  0  rad/m 
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Cylinder  With  Inner  and  Outer  Fluids  of  Air — Young’s  Modulus  Variation.  In  this 
section,  we  compare  the  cylinder/fluid  combination  of  the  previous  section  to  a  cylinder  with  a 
factor-of-10  increase  in  Young’s  modulus.  Figure  103  is  a  simulation  for  a  cylinder  with  a 
Young’s  modulus  of  1  x  10^  Pa  in  contact  with  fluids  of  air,  and  figure  104  is  the  cylinder  (with  a 
Young’s  modulus  of  1  x  10^  Pa)  that  we  have  been  studying  previously. 

The  first  thickness  resonance,  which  was  at  15,198  Hz  (table  11),  has  increased  to  48,330  Hz 
(figure  103).  This  increase  in  frequency  corresponds  to  an  increase  in  the  dilatational  velocity 
for  the  cylinder  (table  12).  The  ratio  of  dilatational  velocities  for  the  two  moduli  is  3.16 
( s  a/IO  ).  The  corresponding  ratio  of  the  first  thickness  resonance  frequencies  is  3.18. 


Table  12.  Young’s  Modulus  and  Dilatational  Phase  Velocity 


Young’s 

Modulus 

(Pa) 

Phase  Velocity  (c/) 

(m/sec) 

1  X  10^ 

447.51 

1  X  10^ 

1,414.9 

The  first  radial  resonance  has  increased  from  1,417  to  4,232  Hz,  as  seen  in  figure  105.  The 
ratio  of  this  frequency  increase  is  2.99,  corresponding  to  a  square  root  of  modulus 
dependence.  From  shell  theory,  a  simple  approximation  for  the  first  radial  resonance  without 
fluid  loading  is 


(326) 


where  a  is  equal  to  the  mean  radius  of  the  cylinder.  This  equation  portrays  the  basic 
relationship  between  the  parameters  for  the  first  radial  resonance.  For  the  nonfluid-loaded  case 
or  the  lightly  fluid-loaded  case  under  consideration  here,  equation  (326)  performs  well  in 
estimating  the  first  radial  resonance  frequency.  Table  13  lists  the  frequencies  calculated  from 
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equation  (326)  and  the  values  taken  from  figure  105. 


Table  13.  Comparison  of  First  Radial  Resonance  Frequency/^ 


Young’s 

First  Radial  Resonance  (Z*^) 

First  Radial  Resonance  if^ 

Modulus 

Equation  (326) 

Figure  105 

(Pa) 

(Hz) 

(Hz) 

lx  10^ 

1,393 

1,417 

1  X  10^ 

r  4,405 

4,232 
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Figure  103.  Radial  Stress  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  WLog{x at  =  1.21  in.  With  F,  =  1  x  10^  Pa,  hy  -  0.6  in., 
n  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  104.  Radial  Stress  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  10Log(T^^^V^i)/^o)^  With  Fj  =  1  x  10*  Pa,  hy  =  0.6  in., 

w  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 
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With  n  =  0,  Diameter  =  3.00  in.,  and  A:  =  0  rad/m 
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Cylinder  With  Inner  and  Outer  Fluids  of  Water — Wall  Thickness  Variation.  In  figures 
106  and  107,  the  pressure  is  calculated  in  the  inner  fluid  at  =  0.01  inch  so  that  the  simulations 
will  be  free  of  fluid  nulls  across  the  frequency  range  under  consideration.  The  thickness 
resonances  are  clearly  visible  in  both  figures  and  correspond  to  cylinder  wall  thicknesses  of  0.3 
and  0.6  inch,  respectively.  Comparison  of  the  responses  is  made  at  A:  =  0  rad/m  in  figure  108. 

The  thickness  resonance  frequencies  are  provided  in  table  14  for  the  simulations  performed 
earlier  with  air  (figures  100  and  101)  and  with  water  (figures  106  and  107).  The  first  thickness 
resonance  is  severely  affected  by  the  presence  of  the  fluid  (water)  in  the  simulation  of  figure 
106.  As  the  cylinder  wall  thickness  is  increased  to  0.6  inch,  the  first  and  third  thickness 
resonance  frequencies  are  not  as  severely  affected;  however,  the  second  thickness  resonance 
frequency  is  reduced  by  50  percent. 


Table  14.  Thickness  Resonance  Table  for  Air  and  Water  With  hi  =  0.3  and  0.6  in. 


Thickness 
Resonance 
Frequency 
Number  (/) 

fti 

Air 

hi  =  0.3  in. 

(Hz) 

fti 

Water 
hi  =  0.3  in. 

(Hz) 

fti 

Air 

hi  =  0.6  in. 

(Hz) 

fti 

Water 
hi  =  0.6  in. 

(Hz) 

1 

30,264 

18,241 

15,198 

13,490 

2 

— 

30,850 

30,532 

20,235 

3 

— 

43,812 

45,221 

41,818 

Cylinder  wall  thickness  variation  also  affects  the  first  radial  resonance,  moving  the  frequency 
from  2816  Hz  at  =  0.3  inch  to  1936  Hz  at  =  0.6  inch.  The  next  section,  which  describes  the 
last  high-frequency  results,  will  examine  a  branch  of  wave  propagation  occurring  near  this 
resonance. 
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Figure  106.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(Pjfri)/F^)2 
at  Vi  =  0.01  in.  With  hi  =  0.3  in.,  n  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  107.  Pressure  Transfer  Surface  Up  to  60  kHz  When  Magnitude  =  10Log(P,(/*i)/P^)2 
at  Vi  -  0.01  in.  With  hi  =  0.6  in.,  n  =  0,  and  Diameter  =  3.00  in.  (Color  Image) 
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Cylinder  With  Inner  and  Outer  Fluids  of  Water — Stress  and  Displacements.  In  this 
last  and  final  section  on  high-frequency  results,  we  shall  examine  the  longitudinal  and  radial 
displacement  fields  that  correspond  to  the  radial  stress  field  at  i\  =  1.21  inches.  For  these 
simulations,  we  return  to  the  original  data  for  the  cylinders  and  fluid  listed  in  tables  2,  4,  and  5, 
where  the  fluids  have  the  properties  of  water  and  the  cylinder  thickness  is  0.3  inch.  The  radial 
stress  is  evaluated  and  displayed  in  figures  109  and  110,  and  the  longitudinal  displacement 
corresponding  to  this  stress  simulation  is  shown  in  figures  111  and  112.  In  a  similar  fashion,  the 
radial  displacement  field  is  displayed  in  figures  113  and  114.  Cuts  through  the  displacement 
surfaces  have  been  made  in  both  wavenumber  and  frequency  and  are  displayed  in  figures  115 
through  118. 

In  figure  115,  the  longitudinal  and  radial  displacements  are  compared  at  2000  Hz.  The 
predominant  displacement  for  the  breathing  wave  is  radial,  and  it  can  be  observed  in  the  figure 
that  the  radial  displacement  is  greater  than  the  longitudinal  displacement  in  the  region  of  the 
breathing  wave  resonance  peak  at  ±170  rad/m.  This  dominance  of  the  radial  displacement  at 
the  breathing  wave  resonance  can  be  observed  in  figures  117  and  118  at  50  and  100  rad/m, 
respectively. 

The  deformation  due  to  the  extensional  wave  is  predominantly  longitudinal.  Therefore,  the 
longitudinal  displacement  will  be  greater  than  the  radial  displacement  in  the  region  of  the 
extensional  wave  resonance  peak.  Figures  117  and  118  illustrate  the  dominance  of  longitudinal 
displacement  at  the  extensional  wave  resonance. 

The  last  feature  in  the  wavenumber-frequency  plane  that  has  not  been  discussed  is  the  branch 
of  wave  propagation  associated  with  the  fluids.  Figure  119  displays  this  region  well;  it  is  a  wire 
frame  representation  of  figure  13,  which  is  a  zoomed-in  view  of  figure  110.  The  phase  velocity 
of  the  dilatational  wave  (also  known  as  the  P-wave  (primary  or  pressure  wave))  for  an  infinite 
media  of  either  inner  or  outer  fluid  is  listed  in  table  1  as  1,500  m/sec.  Close  examination  of  the 
stress  field  indicates  that  there  is  an  increase  in  energy  associated  with  the  fluid  dilatation  wave 
velocity.  The  fluid  branch  is  sharply  resonant  with  very  narrow  extent  in  wavenumber.  The 
cylinder/fluid  combination  considered  here  supports  this  branch  at  the  fluid  P-wave  phase  velocity 
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up  to  the  vicinity  of  the  first  radial  resonance  of  the  cylinder.  At  the  first  radial  resonance  (2816 
Hz),  this  branch  becomes  dispersive,  with  the  phase  and  group  velocities  slowing,  as  shown  in 
table  15.  The  width  in  wavenumber  of  the  resonant  peak  increases  and  a  null  in  the  transmitted 
stress  field  develops  at  the  P-wave  phase  velocity  of  the  fluid  (1500  m/sec).  Above  the  first 
radial  resonance  of  the  cylinder,  the  fluid  P-wave  is  labeled  as  “MODIFIED  FLUID  P-WAVE”  in 
the  figures  that  follow.  The  effect  of  this  fluid/structure  interaction  is  strong  longitudinal 
displacement  in  the  cylinder,  as  evidenced  by  figures  115  through  118. 


Table  15.  Fluid  P-Wave  Phase  and  Group  Velocity  Versus  Frequency  From  Figure  13 


Frequency 

(Hz) 

Phase  Velocity 

(m/sec) 

Group  Velocity  ±75 

(m/sec) 

1,000 

1,500 

1,500 

2,000 

1,500 

1,500 

2,500 

1,500 

1,500 

2,796 

1,460 

1,192 

3,001 

1,426 

1,082 

3,202 

1,416 

1,276 

3,500 

1,393 

1,185 

4,000 

1,296 

872 

4,501 

1,204 

768 

4,990 

1,185 

— 

A  contrast  to  the  previous  case  of  dense  fluids  (water)  is  shown  in  the  simulation  of  figure 
100,  where  the  same  cylinder  was  immersed  in  fluids  of  low  density  (air).  In  this  case,  the  fluid 
P-wave  phase  velocity  is  unaltered  by  the  first  radial  resonance  of  the  cylinder.  The  P-wave 
phase  velocity  of  air  is  343  m/sec  (tables  8  and  9)  and  remains  nondispersive  throughout  the 
wavenumber-frequency  range  of  the  simulation  (figure  1(X)). 
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Figure  109.  Radial  Stress  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  10Log(v^V''i)/^o)^  at  ''i  =  1-21  in.  With  n  ==  0  and 
Diameter  =  3.00  in.  (Color  Image) 


Figure  110.  Radial  Stress  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  lOUogCv^tj^^^yp  )2  ^  ^  21  in.  With  r  =  0  and 

Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  111.  Longitudinal  Displacement  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  10Log(>r^^V^i)/jPo)^  at  =  1-21  in.  With  «  =  0  and  Diameter  =  3.00  in. 

(Color  Image) 


Figure  112.  Longitudinal  Displacement  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  lQiLog{wJ^\r^)IP„Y  at  =  1.21  in.  With  n  =  0  and  Diameter  =  3.00  in. 

(Wire  Frame  Surface) 
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Figure  113.  Radial  Displacement  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  10Log(M^^^(ri)/P^)^  at  =  1.21  in.  With  n  =  0  and  Diameter  =  3.00  in. 

(Color  Image) 


Figure  114.  Radial  Displacement  Transfer  Surface  Up  to  60  kHz  When 
Magnitude  =  10Log(H^‘'^(ri)/F^)2  at  =  1.21  in.  With  n  =  0  and  Diameter  =  3.00  in. 

(Wire  Frame  Surface) 
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Figure  115.  Cylinder  Displacement  Comparison  for  Figures  111  and  113  When 
Magnitude  =  lOLoglX/P^)^  for  Values  of  V  =  ^nd  u^^iry)  at  =  1.21  in. 

With  n  =  0,  Diameter  =  3.00  in.,  and  /=  2000  Hz 
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Figure  116.  Cylinder  Displacement  Comparison  for  Figures  111  and  113  When 
Magnitude  =  lOLoglAT/P^,)^  for  Values  of  V  =  and  u^^(ry)  at  =  1.21  in. 

With  M  =  0,  Diameter  =  3.00  in.,  and  /  =  10000  Hz 


153 


TR  11,067 


Figure  117.  Cylinder  Displacement  Comparison  for  Figures  111  and  113  When 
Magnitude  =  lOLogCZ/P^)^  for  Values  of  V  =  and  at  rj  =  1.21  in. 

With  n  =  0,  Diameter  =  3.00  in.,  and  A:  =  50  rad/m 
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Figure  118.  Cylinder  Displacement  Comparison  for  Figures  111  and  113  When 
Magnitude  =  10Log(X/P(,)^  for  Values  of  Z  =  and  u^^(rx)  at  =  1.21  in. 

With  n  =  0,  Diameter  =  3.00  in.,  and  k  =  100  rad/m 
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Figure  119.  Radial  Stress  Transfer  Surface  With  Magnitude  = 

at  ti  =  1.21  in.  With  n  =  0  and  Diameter  =  3.00  in. 
(Wire  Frame  Surface  of  Figure  13) 
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SINGLE-LAYER  CYLINDER  SIMULATIONS— LONGITUDINAL  SHEAR 
STRESS  EXCITATION 

Up  to  this  point,  the  type  of  excitation  applied  to  the  cylinder  has  been  a  pressure  normal  to 
the  cylinder  surface.  In  this  section,  the  normal  pressure  excitation  will  be  replaced  with  a 
longitudinal  shear  stress  excitation,  once  again  applied  uniformly  over  the  outer  surface  of  the 
cylinder.  The  cylinder  and  fluid  properties  are  taken  from  tables  2,  4  and  5. 

Figures  120  and  121  depict  the  simulation  for  the  n  =  0  excitation.  The  major  branch  of 
wave  propagation  excited  is  the  extensional  wave.  The  phase  velocity  of  the  extensional  wave  is 
roughly  300  m/sec,  and  it  is  nondispersive  over  the  range  of  the  simulation.  The  breathing  wave 
is  faintly  visible  in  figure  121.  In  the  earlier  normal  pressure  simulation  of  figure  20,  we 
observed  a  high-level  breathing  wave  and  a  low-level  extensional  wave;  however,  the 
longitudinal  excitation  creates  a  high-level  extensional  wave  and  a  low-level  breathing  wave. 

The  fluid  P-wave  discussed  in  the  previous  section  is  visible  in  the  n  =  0  simulation.  This 
branch  has  the  same  phase  and  group  velocities  listed  in  table  15. 

The  simulation  for  circumferential  order  number  n  =  1  is  shown  in  figures  122  and  123,  and 
the  simulation  for  «  =  2  is  shown  in  figures  124  and  125.  In  both  cases,  the  second  and  third 
branches  produce  higher  level  responses  than  does  the  first  branch.  The  reason  for  this  is  that 
the  longitudinal  shear  stress  excitation  imparts  energy  more  directly  into  the  longitudinal  and 
circumferential  motions  of  the  cylinder  than  into  the  radial  motion. 

Figure  126  is  a  comparison  of  n  =  0,  1,  and  2  at  50  Hz.  The  major  feature  at  n  =  0  is  the 
extensional  wave,  exhibiting  sharp  roll-off  in  wavenumber  and  a  resonance  amplitude  of 
approximately  20  dB.  At  500  Hz,  in  figure  127,  the  extensional  wave  amplitude  has  dropped  to 
approximately  12  dB  and  the  roll-off  in  wavenumber  is  less  precipitous.  The  fluid  P-wave  is 
indicated  on  figure  128  at  2000  Hz  (the  extensional  wave  of  the  cylinder  is  also  observed  in  this 
region).  In  figure  129  at  5000  Hz,  the  P-wave  has  become  modified  by  the  cylinder,  the  width  of 
the  resonant  peak  has  increased,  and  the  wave  has  become  dispersive.  The  amplitude  of  the 
extensional  wave  has  decreased  in  level. 
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In  each  of  the  comparisons  (figures  126  through  129),  the  n  =  0  excitation  creates  the  largest 
amplitude  pressure  in  the  inner  fluid.  The  n  =  I  and  2  pressure  fields  will  exhibit  strong 
variation  in  amplitude  with  radial  position.  Simulations  of  this  behavior  have  not  been  included 
in  this  report. 
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Figure  120.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/r,)/P^)^  at  r,  =  0.75  in. 

for  n  =  0  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  121.  Pressure  Transfer  Surface  With  Magnitude  =  10Log{Pi(ri)/P^f  at  rj  =  0.75  in. 
for  «  =  0  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  122.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/rj)/P_y)^  at  rj  =  0.75  in. 
for  n  =  1  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  123.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F,(ri)/Pjj.)^  at  =  0.75  in. 
for  n-1  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  124.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P^)^  at  ry  =  0.75  in. 
for  n  =  2  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  125.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F/r,)/P^)2  at  rj  =  0.75  in. 
for  n  =  2  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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SINGLE-LAYER  CYLINDER  SIMULATIONS— CIRCUMFERENTIAL  SHEAR 
STRESS  EXCITATION 

In  this  section,  the  response  of  the  cylinder  to  a  circumferential  shear  stress  excitation  is 
considered.  The  properties  for  the  cylinder  and  fluids  are  the  same  as  those  used  in  the  previous 
section  and  are  taken  from  tables  2, 4  and  5.  Because  there  is  no  circumferential  displacement  in 
an  n  =  0  excitation,  we  will  only  consider  the  n  =  1  and  n  =  2  circumferential  order  numbers  in 
this  section. 

In  figures  130  and  131,  which  depict  the  simulation  for  the  n  =  1  circumferential  shear  stress 
excitation,  the  first  three  branches  of  wave  propagation  are  evident.  The  first  branch  gives  rise  to 
wave  propagation  that  places  the  cylinder  in  the  characteristic  rx  plane-bending  motion.  The 
second  and  third  branches  have  cutoff  frequencies  associated  with  them.  Pressure  is  developed 
in  the  inner  fluid  by  Poisson  coupling  into  a  radial  displacement. 

Figures  132  and  133  display  the  n  =  2  simulation  for  circumferential  shear  stress 
excitation.  Two  branches  of  wave  propagation  are  evident  in  the  simulation,  both  exhibiting 
cutoff  frequencies. 

Figures  134  through  137  are  comparisons  of  the  n  =  1  and  n  =  2  responses  at  various 
frequencies.  For  all  frequencies  and  wavenumbers,  the  n  =  1  excitation  produces  a  larger 
amplitude  pressure  field  in  the  inner  fluid.  The  pressure  field  developed  by  both  the  n  =  1  and 
n  =  2  excitations  is  strongly  dependent  on  radial  position. 
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Figure  132.  Pressure  Transfer  Surface  With  Magnitude  =  lOLog(Pj(r,)/P0)^  at  r,  =  0.75  in. 

for  u  =  2  and  Diameter  =  3.00  in.  (Color  Image) 


Figure  133.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(ri)/Pe)^  at  ri  =  0.75  in. 
for  n  =  2  and  Diameter  =  3.00  in.  (Wire  Frame  Surface) 
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Figure  134.  Comparison  of  Circumferential  Shear  Stress  Excitation  for  Figures  130  and 
132  When  Magnitude  =  lOLogCP/ri^/Pe)^  at  =  0.75  in.  for  Various  n 
With  Diameter  =  3.00  in.  and/=  50  Hz 
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Figure  135.  Comparison  of  Circumferential  Shear  Stress  Excitation  for  Figures  130  and 
132  When  Magnitude  =  lOLog(P/ri)/P0)^  at  =  0.75  in.  for  Various  n 
With  Diameter  =  3.00  in.  and  /  =  500  Hz 
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Figure  136.  Comparison  of  Circumferential  Shear  Stress  Excitation  for  Figures  130  and 
132  When  Magnitude  =  lOLog(F,<ri)/P0)2  at  rj  =  0.75  in.  for  Various  n 
With  Diameter  =  3.00  in.  and/=  2000  Hz 
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Figure  137.  Comparison  of  Circumferential  Shear  Stress  Excitation  for  Figures  130  and 
132  When  Magnitude  =  lOLogCP/rij/Pe)^  at  =  0.75  in.  for  Various  n 
With  Diameter  =  3.00  in.  and  /=  5000  Hz 
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TWO-LAYER  CYLINDER  SIMULATIONS 

Thus  far  we  have  considered  simulations  of  a  cylinder  composed  of  a  single  homogeneous 
layer  of  material  in  contact  with  fluids  on  the  inner  and  outer  surfaces.  In  this  section,  we  will 
employ  the  full  model  developed  in  this  report  for  the  two-layer  cylinder  with  fluids.  The 
system  matrix  is  given  by  equation  (150),  with  elements  given  by  equations  (151)  through 
(319).  Material  properties  and  dimensions  for  the  two-layer  cylinder  are  shown  in  table 
16.  The  first  or  inner  cylinder  has  material  properties  similar  to  hard  plastic,  such  as 
polycarbonate.  The  outer  cylinder  (cylinder  2)  has  material  properties  similar  to  those  of 
urethane  or  the  rubber  used  for  automobile  tires.  The  inner  fluid  properties  are  consistent  with 
air  and  are  listed  in  table  17.  The  outer  fluid  material  properties  are  those  of  water  and  are  given 
in  table  18. 


Table  16.  Two-Layer-Cylinder  Properties — 0.650-in.  Diameter 


Cylinder  #1  Property 

Cylinder  #2  Property 

Definition 

9  N 

E.  =  2.4x10  ^ 

7  N 

=  1.5x10  ^ 

^  m2 

Young’s  Modulus 

11 

P 

^2  =  0.20 

Structural  Loss  Factor 

Pi  =  2113^ 

P2  =  1000^ 

Density 

Vj  =  0.40 

V2  =  0.45 

Poisson’s  Ratio 

a  =  0.160in. 

b  =  0.260in. 

Inner  Radius 

b  =  0.260in. 

c  =  0.325in. 

Outer  Radius 

=  O.lOOin. 

h2  =  0.065  in. 

Thickness 

The  model  under  consideration  is  composed  of  four  distinct  layers.  The  simulations  will 
contain  contributions  of  wave  propagation  and  interaction  from  each  of  the  four  layers.  The 
resulting  composite  system  will  produce  branches  of  wave  propagation  that  are  related  to  the 
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dilatational  and  transverse  waves  of  the  media  but  do  not  necessarily  propagate  at  the  same 
velocity.  The  velocities  listed  in  table  19  will  serve  as  a  guide  for  analyzing  the  responses  that 
follow.  The  simulations  for  the  two-layer  cylinder  will  be  displayed  as  color  images  using  the 
red/yellow/green/blue  color  map  as  opposed  to  the  red  temperature  color  map  that  was  previously 
used  for  the  single-layer  cylinder. 


Table  17.  Two-Layer-Cylinder  Inner  Fluid  Properties  (Air) 


Property 

Definition 

P,  =  1-20^ 

Density 

c.  =  343.0  ^ 

'  sec 

Velocity  of  Sound 

Table  18.  Two-Layer-Cylinder  Outer  Fluid  Properties  (Water) 


Property 

Definition 

p.  =  1000.0^ 

Density 

c  =  1500.0  ^ 
sec 

Velocity  of  Sound 

Table  19.  Two-Layer-Cylinder/Fluid  Phase  Velocities 


Inner  Fluid 
Phase  Velocity 

(m/sec) 

Cylinder  #1 
Phase  Velocity 

(m/sec) 

Cylinder  #2 
Phase  Velocity 

(m/sec) 

Outer  Fluid 
Phase  Velocity 

(m/sec) 

Cl 

— 

1,560.0 

238.5 

— 

Ct 

— 

636.9 

71.9 

— 

Ce 

— 

1,065.8 

122.5 

— 

Ci 

343.0 

— 

— 

— 

— 

— 

— 

1,500.0 
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Radial  Pressure  Excitation 

The  first  series  of  simulations  are  for  a  radial  pressure  excitation  applied  to  the  outer  surface 
of  the  cylinder,  as  shown  in  figure  1 .  The  circumferential  strain  is  the  quantity  calculated  just 
below  the  outer  surface  at  =  0.324  inch,  and  it  is  displayed  as  a  color  image  using  the  red/ 
yellow/green/blue  color  map  in  figure  138  and  as  a  wire  frame  surface  in  figure 
139.  Comparison  of  the  amplitude  in  frequency  is  made  in  figure  140. 

The  first  radial  resonance  frequency  occurs  at  6210  Hz  for  this  system.  The  P-wave  in  the 
inner  fluid  is  faintly  visible  as  a  series  of  dots  at  a  phase  velocity  of  343  m/sec  in  figure 
138.  This  wave  is  unaffected  by  the  structure,  which  is  consistent  with  the  previous  observation 
of  what  occurred  when  air  was  used  as  the  fluid  with  the  single-layer  cylinder.  The  P-wave  of 
the  outer  fluid  experiences  dispersion  due  to  interaction  with  the  cylinder  above  the  radial 
resonance  frequency.  The  other  major  branch  of  wave  propagation  evident  in  the  simulation  is 
an  extensional  wave  for  the  composite  cylinder.  The  phase  velocity  of  this  wave  is  shown  in 
table  20;  its  value  falls  between  those  listed  in  table  19  for  the  individual  layers. 


Table  20.  Extensional  Wave  Phase  Velocity  From  Figure  138 


Frequency 

(Hz) 

Phase  Velocity 

(m/sec) 

1,000 

800 

2,000 

800 

3,000 

800 

4,000 

800 

5,000 

800 

6,000 

800 

7,000 

700 

7,500 

663 

7,800 

610 
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Figures  141  through  143  are  simulations  of  the  circumferential  strain  at  =  0.259  inch,  which 
is  just  below  the  outer  radius  of  the  inner  cylinder.  The  effect  of  the  inner  fluid  P-wave  on  the 
strain  field  is  captured  in  the  6000-Hz  curve  in  figure  143  (note  the  small  bump  at  approximately 
100  rad/m).  This  P-wave  produces  a  very  narrow  spike  in  wavenumber.  The  wavenumber 
interval  chosen  for  the  simulations  is  not  small  enough  to  capture  the  P-wave  continuously,  which 
is  why  it  appears  as  a  series  of  dots  in  the  color  images.  Decreasing  the  wavenumber  interval 
would  lead  to  a  corresponding  increase  in  calculation  time  and  was  not  considered  important 
enough  for  the  purposes  of  this  report. 

The  circumferential  strain  field  is  evaluated  at  =  0.250  inch  (which  is  within  the  confines  of 
the  thickness  of  the  inner  cylinder),  as  shown  in  figures  144  through  146.  In  figures  147  through 
149,  the  circumferential  strain  field  is  evaluated  at  rj  =  0.161  inch,  which  is  just  inside  the  inner 
radius  of  the  inner  cylinder.  Comparison  at  ^  =  0  is  made  in  figure  150  for  the  four  radial 
locations  previously  mentioned,  and  the  radial  resonance  can  be  clearly  observed  in  all  the 
curves.  Figures  151  through  154  are  comparisons  at  various  frequencies  for  each  of  the  four 
radial  locations.  Although  the  extensional  wave  in  the  cylinder  and  the  P-wave  in  the  outer  fluid 
have  been  labeled  only  once  in  figures  153  and  154  (at  =  0.32  inch)  to  prevent  cluttering  of  the 
plots,  they  can  also  be  seen  at  the  radial  locations  shown  on  the  other  curves.  The  P-wave  in  the 
inner  fluid  is  evident  in  figure  153.  This  wave  produces  a  narrow  spike  in  wavenumber,  which 
was  not  captured  by  the  calculation  in  the  negative  wavenumber  region  because  the  algorithms  do 
not  necessarily  evaluate  the  equations  at  precisely  the  same  wavenumber  on  both  sides  of  zero. 

The  attenuation  of  circumferential  strain  sensitivity  is  displayed  in  figure  155.  The 
magnitude  plotted  in  this  figure  is  the  difference  in  level  between  rj  =  0.259  inch  and  =  0.324 
inch.  The  curves  indicate  the  wavenumber-filtering  property  of  the  structure  due  to  normal 
pressure  excitation  from  just  beneath  the  outer  radius  of  the  outer  cylinder  to  just  beneath  the 
outer  radius  of  the  inner  cylinder. 
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Figure  138,  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  lOLog(800^^(r,)/P„)^  at  Tj  =  0.324  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Color  Image) 


Figure  139.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  10Log(£ee<^2(ri)/P^)2  at  rj  =  0.324  in.  With  «  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  141.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  lOLoglege^VuV^o)^  at  ri  =  0.259  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Color  Image) 


Figure  142.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  lOLogCegg^Vu)/^,,)^  at  G  =  0.259  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  144.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  at  =  0.250  in.  With  «  =  0  and 

Diameter  =  0.650  in.  (Color  Image) 


Figure  145.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  lOLog(e0e‘^VGU^o)^  at  Vy  =  0.250  in.  With  «  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 


181 


Figure  148.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  lOLog(8e0^^(ri;/Po)2  at  =  0.161  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  150.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  138, 141, 144, 
and  147  When  Magnitude  =  10Log(89e^Y^i>/Fo)^  at  Various  ti  With  /  =  1  or  2,  n  =  0, 

Diameter  =  0.650  in.,  and  k  =  0  rad/m 
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Figure  151.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  138, 141, 144, 
and  147  When  Magnitude  =  lOLogleee^^Vil/Ro)^  at  Various  With  i  =  1  or  2,  n  =  0, 

Diameter  =  0.650  in.,  and  /  =  1000  Hz 
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Figure  153.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  138, 141, 144, 
and  147  When  Magnitude  =  lOLog(ee0^Y^ij/jPo)^  at  Various  With  i  =  1  or  2,  n  =  0, 

Diameter  =  0.650  in.,  and  /=  6000  Hz 
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Figure  155.  Attenuation  of  Circumferential  Strain  From  Figures  138  and  141 
When  Magnitude  =  (lOLog(E0e^Vo.259)/Pj2 .  iOLog(80e^^fO.324)/P„)2)  at  Various/ 

With  «  =  0  and  Diameter  =  0.650  in. 
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Longitudinal  and  Radial  Displacements.  The  longitudinal  and  radial  displacements  that 
correspond  to  the  location  =  0.259  inch  (as  discussed  in  the  previous  section)  are  now 
calculated.  The  simulation  for  longitudinal  displacement  is  shown  in  figures  156  and  157;  the 
simulation  for  radial  displacement  is  presented  in  figures  158  and  159.  Comparisons  of  these 
two  surfaces  are  made  at  selected  frequencies  in  figures  160  through  163. 

At  1000  Hz  in  figure  160  and  3000  Hz  in  figure  161,  the  extensional  wave  is  the  major  branch 
of  wave  propagation.  Moving  up  to  6000  Hz  in  figure  162  and  7500  Hz  in  figure  163,  the  inner 
and  outer  fluid  P-waves  become  evident.  Because  the  wavenumber  increment  used  for  the 
simulation  is  larger  than  the  width  in  wavenumber  of  the  inner  fluid  P-wave  resonance,  the  inner 
fluid  P-wave  shows  up  intermittently  in  the  simulations. 
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Figure  156.  Longitudinal  Displacement  Transfer  Surface  With 
Magnitude  =  10hog(w at  =  0.259  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Color  Image) 


Figure  157.  Longitudinal  Displacement  Transfer  Surface  With 
Magnitude  =  10Log{w J^^(ri)/Pg)^  at  =  0.259  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  158.  Radial  Displacement  Transfer  Surface  With  Magnitude  =  ySLog{Uc^\rx)IP 
at  =  0.259  in.  With  n  =  0  and  Diameter  =  0.650  in.  (Color  Image) 


Figure  159.  Radial  Displacement  Transfer  Surface  With  Magnitude  =  10Log(Mj.^'(rD/P(,) 
at  =  0.259  in.  With  n  =  0  and  Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  160.  Comparison  of  Displacement  Distribution  for  Figures  156  and  158  When 
Magnitude  =  10Log(X/Po)2  for  Values  of  X  =  and  at  =  0.259  in. 

With  n  =  0,  Diameter  =  0.650  in.,  and/=  1000  Hz 
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Figure  161.  Comparison  of  Displacement  Distribution  for  Figures  156  and  158  When 
Magnitude  =  10Log(Z/P^)2  for  Values  of  V  =  and  at  =  0.259  in. 

With  /I  =  0,  Diameter  =  0.650  in.,  and  /  =  3000  Hz 
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Figure  162.  Comparison  of  Displacement  Distribution  for  Figures  156  and  158  When 
Magnitude  =  10Log(A/P<,)^  for  Values  of  V  =  and  u^\r^  at  =  0.259  in. 

With  «  =  0,  Diameter  =  0.650  in.,  and/=  6000  Hz 
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Figure  163.  Comparison  of  Displacement  Distribution  for  Figures  156  and  158  When 
Magnitude  =  10Log(X/P^)2  for  Values  of  Z  =  and  at  =  0.259  in. 

With  /I  =  0,  Diameter  =  0.650  in.,  and/=  7500  Hz 
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Longitudinal  Shear  Stress  Excitation 

In  this  section,  the  two-layer  cylinder  is  excited  with  a  longitudinal  shear  stress 
excitation.  The  fluids  are  the  same  as  those  defined  previously  in  tables  17  and  18.  The 
circumferential  strain  sensitivity  is  evaluated  at  =  0.324  inch  in  figures  164  and  165.  This 
radial  location  is  just  below  the  outer  surface  of  the  second  cylinder.  Cuts  through  the  surface 
are  made  in  figure  166  at  various  frequencies.  The  major  branch  of  wave  propagation  in  this 
simulation  for  the  composite  cylinder  is  the  extensional  wave  whose  amplitude  decreases  as 
frequency  increases  (see  figure  166). 

At  a  depth  of  rj  =  0.259  inch,  just  below  the  outer  radius  of  the  inner  cylinder,  the  simulation 
for  circumferential  strain  is  presented  in  figures  167  and  168.  Cuts  through  this  surface  at 
various  frequencies  are  presented  in  figure  169.  As  was  the  case  at  rj  =  0.324  inch,  the 
extensional  wave  is  the  dominant  branch  of  wave  propagation,  although  the  decrease  in  its 
amplitude  with  frequency  is  greater  at  rj  =  0.259  inch  than  at  rj  =  0.324  inch. 

Comparisons  of  the  circumferential  strain  at  the  two  radial  locations  are  made  in  figures  170 
through  174.  Above  the  extensional  wavenumber,  there  is  significant  attenuation  of 
circumferential  strain  sensitivity.  At  +300  rad/m,  this  sensitivity  has  been  reduced  by 
approximately  30  dB.  The  extensional  wave  resonance  peak  level,  which  decreases  as  rj 
decreases  from  0.324  to  0.259  inch,  can  be  observed  in  figures  171  through  174. 

The  attenuation  of  circumferential  strain  sensitivity  is  displayed  in  figure  175.  The 
magnitude  plotted  in  this  figure  shows  the  difference  in  level  between  =  0.259  inch  and 
rj  =  0.324  inch.  The  curves  reveal  the  wavenumber-filtering  property  of  the  structure  due  to 
longitudinal  shear  stress  excitation  from  just  beneath  the  outer  radius  of  the  outer  cylinder  to  just 
beneath  the  outer  radius  of  the  inner  cylinder. 

A  comparison  of  figure  175  with  figure  155  reveals  the  different  amounts  of  filtering  that  the 
structure  presents  to  the  two  forms  of  excitation  (i.e.,  longitudinal  shear  stress  and  normal 
pressure  excitation). 
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Figure  164.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  10Log(89e^^(ri)/P^)2  at  ri  =  0.324  in.  With  /i  =  0  and 
Diameter  =  0.650  in.  (Color  Image) 


Figure  165.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  10Log(8ee^^(ri)/P;^)2  at  rj  =  0.324  in.  With  «  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  167.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  lOLogCeee^^frD/P^^)^  at  ri  =  0.259  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Color  Image) 


Figure  168.  Circumferential  Strain  Transfer  Surface  With 
Magnitude  =  10Log(ee9^^(ri;/PD^  at  =  0.259  in.  With  n  =  0  and 
Diameter  =  0.650  in.  (Wire  Frame  Surface) 
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Figure  170.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  164  and  167 
When  Magnitude  =  lOLog(e00^Y^t)/F;c)^  at  Various  With  i  =  1  or  2,  n  =  0, 
Diameter  =  0.650  in.,  and  k  =  0  rad/m 


Figure  172.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  164  and  167 
When  Magnitude  =  lOLog(8e0^Vt)/F;c)^  at  Various  With  i  =  1  or  2,  n  =  0, 
Diameter  =  0.650  in.,  and/=  3000  Hz 
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Figure  173.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  164  and  167 
When  Magnitude  =  lOLog(ee0^Vi)/P;f)^  at  Various  With  i  =  1  or  2, «  =  0, 
Diameter  =  0.650  in.,  and/=  6000  Hz 
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Figure  175.  Attenuation  of  Circumferential  Strain  From  Figures  164  and  167 
When  Magnitude  =  (lOLog(ee0^^fO.259)/F_^)2  -  10Log(eee^^(0.324)/P_^)2)  at  Various/ 

With  n  =  0  and  Diameter  =  0.650  in. 


m 


TR  11,067 


Comparison  of  Radial  Pressure  and  Longitudinal  Shear  Stress  Excitations.  Previously, 
we  have  examined  the  circumferential  strain  sensitivity  to  normal  pressure,  P^,  and  to 
longitudinal  shear  stress,  In  this  section,  we  will  compare  this  sensitivity  to  P^  (figure  141) 
and  P^  (figure  167)  at  various  frequencies.  Figures  176  through  181  are  cuts  in  frequency  for 
this  comparison. 

At  40  Hz  (figure  176),  the  extensional  wave,  at  low  wavenumber,  shows  the  circumferential 
strain  sensitivity  to  longitudinal  shear  stress  to  be  above  the  sensitivity  to  normal  pressure  by 
40  dB.  As  frequency  increases,  the  magnitude  of  the  extensional  wave  generated  by  the 
longitudinal  shear  stress  excitation  decreases,  and  at  6000  Hz  the  difference  between  the  two 
responses  in  the  region  of  the  extensional  wave  is  approximately  13  dB. 

In  general,  the  figures  show  that  the  circumferential  strain  sensitivity  to  longitudinal  shear 
stress  is  lower  than  the  normal  pressure  sensitivity  above  100  rad/m.  The  circumferential  strain 
sensitivity  to  normal  pressure,  unlike  longitudinal  shear  stress,  is  relatively  flat,  with  increasing 
wavenumber  above  the  cylinders  extensional  wave  region. 
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Figure  177.  Comparison  of  Circumferential  Strain  Distribution  for  Figures  141  and  167 
When  Magnitude  =  10Log(eee^V^i)/Fj)2  With  i  =  o  or  x,  n  =  0, 

Diameter  =  0.650  in.,  and  /  =  500  Hz 
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CONCLUSIONS  AND  RECOMMENDATIONS 

An  exact  closed-form  solution  to  the  problem  of  forced  harmonic  excitation  of  a  two-layer 
cylinder  with  inner  and  outer  fluids  has  been  derived.  Results  are  presented  for  the  dynamic 
response  of  a  single-layer  cylinder  with  inner  and  outer  fluids  of  water  and  for  a  two-layer 
cylinder  with  air  as  the  inner  fluid  and  water  as  the  outer  fluid. 


Single-Layer  Summary 

The  elasticity  model  is  compared  to  membrane  and  bending  shell  models  developed 
previously.  Comparisons  among  the  models  are  made  for  a  3.00-inch-diameter  cylinder  and  a 
0.670-inch-diameter  cylinder.  To  summarize  the  results,  for  the  3.00-inch-diameter  cylinder,  the 
shell  models  are  in  good  agreement  with  the  elasticity  model  for  the  region  extending  to  1000  Hz 
and  ±75  rad/m.  Outside  of  this  region,  the  difference  increases  beyond  5  dB,  and  the  elasticity 
solution  should  be  used.  For  the  n  =  0  case,  the  fluid  loading  influences  the  first  radial  natural 
frequency  of  the  shell  models  so  dramatically  that  the  error  (relative  to  the  elasticity  solution)  is 
greater  than  a  factor  of  three  as  evidenced  by  figure  22.  For  the  n  =  2  case,  the  membrane  model 
should  not  be  used  at  all  beeause  it  misrepresents  the  behavior  of  the  first  branch  (i.e.,  lack  of 
cutoff  frequency  and  appropriate  branch  structure).  Although  not  shown  in  this  report,  similar 
behavior  will  occur  with  the  membrane  shell  model  for  «  >  2 . 

For  the  0.670-inch-diameter  cylinder,  the  membrane  shell  models  are  generally  within  4  dB  of 
the  elasticity  model  for  the  n  =  0  circumferential  order  number  in  the  region  ±150  rad/m  and  0  to 
5000  Hz.  If  it  is  necessary  to  obtain  a  simulation  with  an  error  of  less  than  4  dB,  the  elasticity 
model  should  be  used.  The  shell  models  should  not  be  used  for  the  determination  of  the  first 
radial  natural  frequency  of  this  cylinder/fluid  combination;  the  elasticity  model  should  be 
used.  For  the  higher  circumferential  order  numbers,  it  should  be  remembered  that  although  the 
membrane  model  will  not  support  a  cutoff  frequency  for  the  first  branch  of  wave  propagation,  the 
bending  shell  model  will  be  within  4  dB  of  the  elasticity  solution.  If  an  error  of  less  than  4  dB  is 
required,  the  elasticity  model  should  be  used. 
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The  pressure  field  in  the  inner  fluid  consists  of  propagating  and  nonpropagating 
regions.  Below  the  wavenumber  w/c,-,  the  pressure  field  propagates  and  is  given  by  the  Bessel 
function  of  equation  (67).  Above  w/c,  the  pressure  field  becomes  increasingly  evanescent,  with 
increasing  wavenumber  and  decreasing  radial  position  r,  and  is  given  by  equation  (71).  This 
behavior  serves  to  attenuate  breathing  wave  and  extensional  wave  pressure  levels  (which  are 
locally  generated  at  the  cylinder/fluid  surface)  as  the  measurement  point  moves  toward  the  center 
of  the  cylinder  (r  =  0)  or  to  a  higher  wavenumber. 

Decreasing  the  density  of  the  inner  fluid  lowers  the  first  radial  resonance  of  the  cylinder  fluid 
combination.  Decreasing  the  density  also  increases  the  phase  velocity  of  the  breathing 
wave.  A  factor  of  four  decrease  in  inner  fluid  density  almost  doubles  the  breathing  wave  phase 
velocity  from  66  to  130  m/sec. 

A  comparison  of  shell  and  elasticity  models  at  high  frequencies  (up  to  60  kHz)  shows  that  the 
shell  models  do  not  accurately  resolve  the  behavior  of  the  cylinder.  For  the  cases  considered 
here  (which  represent  heavy  fluid  loading),  the  elasticity  model  accurately  predicts  the  first  radial 
natural  frequency  of  the  system  (the  shell  models  do  not).  Furthermore,  in  contrast  to  the  shell 
model,  the  elasticity  model  reveals  the  specific  dynamics  associated  with  the  thickness  of  the 
cylinder.  Dilatation  wave  propagation  creates  standing  waves  across  the  thickness  of  the 
cylinder  at  frequencies  (/^,)  corresponding  to  an  integer  number  of  half  wavelengths  for  light  fluid 
loading  (air).  Doubling  the  cylinder  wall  thickness  halves  the  frequencies  (tables  10  and  11), 
and  increasing  the  Young’s  modulus  increases  the  resonant  frequencies,  as  expected  (table 
13).  For  heavy  fluid  loading  (water),  the  frequencies  that  result  deviate  from  the  integer- 
number-of-half-wavelength  relationship  observed  when  air  is  used  as  the  inner  and  outer  fluids 
(table  14).  Nulls  occur  in  the  inner  fluid  pressure  field  due  to  the  nature  of  the  field  (the  null 
frequency  is  given  by  equation  (325)  and  is  a  function  of  the  dilatational  phase  velocity  of  the 
inner  fluid,  the  radial  position  rj,  and  the  wavenumber  k).  The  first  two  null  frequencies  are 
tabulated  in  table  7  for  the  3.00-inch  cylinder  with  water  as  the  internal  fluid. 

The  dilatational  or  P-wave  phase  velocity  of  the  inner  and  outer  fluids  interacts  with  the 
cylinder.  For  the  3.00-inch-diameter  cylinder  with  fluids  of  water  (figures  109  and  119),  close 


218 


TR  11,067 


examination  of  the  stress  field  indicates  that  there  is  an  increase  in  energy  associated  with  the  fluid 
dilatation  wave  velocity;  that  is,  the  fluid  branch  becomes  sharply  resonant  with  a  very  narrow 
extent  in  wavenumber.  The  combination  of  cylinder  and  fluids  considered  here  supports  this 
branch  at  the  fluid  P-wave  phase  velocity  up  to  the  vicinity  of  the  first  radial  resonance  of  the 
cylinder.  At  the  radial  resonance  (2816  Hz),  this  branch  becomes  dispersive,  with  the  phase  and 
group  velocities  slowing,  as  shown  in  table  15.  The  width  in  wavenumber  of  the  resonant  peak 
increases  and  a  null  in  the  transmitted  stress  field  develops  at  the  P-wave  phase  velocity  of  the 
fluid  (1,500  m/sec).  Above  the  first  radial  resonance  of  the  cylinder,  the  fluid  P-wave  is 
modified  by  the  presence  of  the  cylinder  and  becomes  dispersive. 

In  contrast  to  the  previously  discussed  case  of  dense  fluids  (water)  is  the  simulation  shown  in 
figure  100,  where  the  3.00-inch  cylinder  was  immersed  in  fluids  of  low  density  (air).  In  this 
case,  the  fluid  P-wave  phase  velocity  is  unaltered  by  the  first  radial  resonance  of  the 
cylinder.  The  P-wave  phase  velocity  of  air  is  343  m/sec  (tables  8  and  9)  and  remains 
nondispersive  throughout  the  wavenumber-frequency  range  of  the  simulation  (figure  100). 

Normal  pressure  excitation  {n  =  0)  applied  to  the  cylinder  creates  a  high-amplitude  breathing 
wave  and  a  low-amplitude  extensional  wave;  however,  the  longitudinal  shear  stress  excitation 
creates  a  high-amplitude  extensional  wave  and  a  low-amplitude  breathing  wave.  The  amplitude 
of  the  extensional  wave  generated  by  the  longitudinal  shear  stress  excitation  decreases  with 
increasing  frequency.  The  interaction  of  the  fluid  P-wave  and  the  cylinder  during  longitudinal 
shear  stress  excitation  («  =  0)  parallels  that  of  the  radial  pressure  excitation  (i.e.,  nondispersive 
behavior  up  to  the  first  radial  resonance  frequency  of  the  cylinder  and  dispersive  behavior  above 
the  resonance  frequency). 

The  n  =  1  circumferential  order  number  excitation  (normal  pressure,  longitudinal  shear  stress, 
and  circumferential  shear  stress)  generates  a  branch  of  wave  propagation  that  causes  the  cylinder 
to  undergo  an  rx  plane-bending  deformation,  as  depicted  in  figure  9.  Cutoff  frequencies  are 
associated  with  the  second  and  third  branches  of  wave  propagation  for  n  =  1  and  n  =  2  for  both 
longitudinal  and  circumferential  shear  stress  excitations. 
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Two-Layer  Summary 

The  two-layer  cylinder  simulations  (circumferential  strain  with  normal  pressure  excitation) 
are  performed  with  air  as  the  inner  fluid  and  water  as  the  outer  fluid.  The  first  radial  resonance 
frequency  for  the  composite  cylinder  and  fluids  occurs  at  6210  Hz.  The  P-wave  of  the  inner 
fluid  is  nondispersive  throughout  the  range  of  the  simulations.  The  P-wave  of  the  outer  fluid  is 
nondispersive  up  to  the  first  radial  resonance  of  the  cylinder  and  then  becomes  dispersive.  The 
extensional  wave  phase  velocity  for  the  two-layer  cylinder  is  800  m/sec,  becoming  dispersive 
above  the  cylinder  radial  resonance  frequency  of  6210  Hz.  Figure  155  displays  the  attenuation 
of  circumferential  strain  across  the  outer  layer  of  the  composite  cylinder. 

A  circumferential  strain  excitation  (n  =  0)  applied  to  the  composite  cylinder  produces  an 
extensional  wave  whose  amplitude  is  20  dB  larger  at  1000  Hz  and  5  dB  larger  at  7500  Hz  than  the 
extensional  wave  generated  by  the  normal  pressure  excitation.  The  outer  fluid  P-wave  is  evident 
in  the  simulations  using  the  circumferential  shear  stress  excitation  at  approximately  the  same 
amplitude  as  for  the  normal  pressure  excitation. 

These  results  provide  a  small  sampling  of  the  problems  that  can  be  analyzed  with  the 
mathematical  models  developed  in  this  report.  It  should  be  emphasized  that  the  elasticity 
models  represent  the  “truth”  for  the  structures  (subjected  to  harmonic  forcing  functions  in  space 
and  time)  presented  here.  There  are  no  simplifying  assumptions  that  discard  information  as  in 
shell  analysis.  The  elasticity  models  can  be  used  for  analyzing  structures  that  range  from  the 
dimensions  of  a  human  hair  to  the  largest  cylinder  that  would  be  practical  to  manufacture  or 
would  occur  in  nature.  Furthermore,  the  frequency  range  allowed  by  these  models  is 
exceedingly  broad.  In  other  words,  the  simulations  presented  here  could  be  said  to  explore  only 
the  lower  limit  of  what  is  possible  to  study  with  the  elasticity  models.  The  primary  challenge 
becomes  a  numerical  one,  where  the  limitations  occur  in  the  ability  to  compute  the  Bessel 
functions,  invert  the  system  matrix,  and  maintain  enough  precision  throughout  the  entire 
calculation  to  ensure  a  valid  solution  to  the  equations. 
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APPENDIX  A 

BESSEL  FUNCTIONS  OF  COMPLEX  ARGUMENT,  SERIES 

REPRESENTATION 

The  inclusion  of  structural  damping  in  the  cylinder  necessitates  the  evaluation  of  Bessel 
functions  of  complex  argument.  The  series  representation  of  these  functions  is  then  required. 
From  Abramowitz  and  Stegun,^  equations  (A-1)  through  (A-6)  are  used  for  complex  argument  z, 
where  z  ranges  from  0  <  z  <  3 .  The  circumferential  order  number  n  appears  in  the  equations  as 
well. 

The  Gamma  function  is  given  by 

r(n  +  .fi:+l)  =  {n+K)\.  (A-1) 


o 

The  Psi  or  Digamma  function  is 

n  -  1 

\(r(l)=-y,  xi/(n)  = -7+  ^  A:-i,  n>2  . 

k=  1 

The  value  used  for  Euler’s  constant  is  7=  0.577215664901532860606512. 


The  Bessel  function  of  the  first^  kind  is 


(A-2) 


CO 


k  =  0 


/rir  (/I  +  k 


(A-3) 


1.  M.  Abramowitz  and  1.  A.  Steguriy  Handbook  of  Mathematical  Functions,  U.S.  Government  Printing  Office, 
Washington,  DC,  June  1964. 

2.  M.  Abramowitz  and  1.  A.  Stegun,  p.  255,  equation  (6.1.6). 

3.  M.  Abramowitz  and  I.  A.  Stegun,  p.  258,  equation  (6.3.2). 

4.  M.  Abramowitz  and  I.  A.  Stegun,  p.  360,  equation  (9.1.10). 
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The  Bessel  function  of  the  second^  kind  over  the  range  0  <  z  <  3  is 


111 

f  n-\ 

{h  ^  k  —  1 ) ! 

71 

V 

y 

^  k\ 

^k  =  0 

A 

/ 


k  2, 

+  -In 

n 


k  =  0 


(\|/ (/:+!)  +  \|/(n  +  A:+l)) 


/cl(n  +  k)r 


(A-4) 


The  Hankel  function^  of  the  first  kind  is  given  by 


(A-5) 


The  expression^  used  for  the  modified  Bessel  function  /„  over  the  range  0  <  z  <  3  is 


4(2) 


/c  =  0 


k!r(n  +  k+  1)  ■ 


(A-6) 


The  expression^  used  for  the  modified  Bessel  function  over  the  range  0  <  z  <  3  is 


if 

2V2 


^  -  1 

Z 

Vjt  =  0 


(n-k  -  1)1 

k! 


-iy)'+(-l)”"nn(y/„(z) 


Z 

+  (-1)1 


^  (\\f(k+l) +\^in  +  k+l)) 


X  I  1  2^'^ 

A  l-z2 


k\  {n  +  k)\ 


(A-7) 


5.  M.  Abramowitz  and  I.  A.  Stegun,  p.  360,  equation  (9.1.11). 

6.  M.  Abramowitz  and  I.  A.  Stegun,  p.  358,  equation  (9.1.3). 

7.  M.  Abramowitz  and  I.  A.  Stegun,  p.  375,  equation  (9.6.10). 

8.  M.  Abramowitz  and  I.  A.  Stegun,  p.  375,  equation  (9.6.1 1). 
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For  the  argument  range  3  <  z  <  oo ,  asymptotic  expansions  are  used  from  Kom  and  Kom:^ 


= 


A„  (z)  cos  z 
TIZV  "  V  2 


nit  7t  (  nn  n 


(A-8) 


and 


(A-9) 


where 


4  (4^2- 1)  (4^2 -9)  I  (4n2_  1)  (4n2_9)  (4^2_25)  (4^2_49) 

”  2!(8z)2  4!(8z)4 


and 


4n^  -  1  (4n^  -  1)  -  9)  {4n'^  -  25)  ^ 

8z  3!(8z)3 


(A-10) 


9.  G.  A.  Kom  and  T.  M.  Kom,  Mathematical  Handbook  for  Scientists  and  Engineers,  McGraw-Hill  Book 
Company,  New  York,  pp.  868-869,  equations  (21.8-44)  and  (21.8.45). 
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APPENDIX  B 

DERIVATIVES  OF  BESSEL  FUNCTIONS 

An  expression  for  the  first  derivative  of  the  Bessel  functions  needed  for  the  evaluation  of  the 
displacement  potentials  is  taken  from  Arpaci,^  and  is  seen  here  in  equations  (B-1)  and  (B-2): 


^iZ^(ax))  = -aZ„^i(ax) +^Z„(ax),  Z  = 

^  (B-1) 

and 

^ (/„  (ax) )  =  a/„ ^  1  (ax)  -t-  (ax)  .  (B-2) 

The  argument  of  the  Bessel  functions  is  composed  of  a  (a  constant)  and  x  (a  complex  number). 

An  expression  for  the  second  derivative  of  the  Bessel  functions  listed  in  equation  (B-1)  is 
derived  by  differentiating  this  equation  again  with  respect  to  x,  which  requires  the  use  of  the  prod¬ 
uct  rule  on  the  second  term  of  the  right-hand  side  of  equation  (B-1)  and  the  following  effect  of  the 
recursion  relationships  for  Bessel  functions: 

=  -aZ^^^iax)  +^z^^^(ax),  Z  = 

^  .  (B-3) 


The  final  expression  for  the  second  derivative  is 

^ZJax)  =  ^(Z„_2(ax)  -  (2Z„(ax)  -(-Z„^2 ) ) 

Z  =  j,y,k,h^^\h^^^ 


(B-4) 


1.  V.  S.  Arpaci,  Conduction  Heat  Transfer,  Addison-Wesley  Publishing  Company,  Reading,  MA,  p.  139, 
equation  (3-139). 


B-l/B-2 
Reverse  Blank 


TR  11,067 


APPENDIX  C 

OUTPUT  QUANTITY  EQUATIONS 

The  equations  necessary  to  compute  the  quantities  chosen  for  the  simulations  are  listed  in  this 
appendix.  The  system  matrices,  equations  (85)  and  (150),  are  inverted  for  one  excitation  at  a 
time  (i.e.,  P^,  or  Pq).  The  fluid  displacement  potential  constant  D  is  used  with  equations  (67) 

and  (71)  to  compute  the  magnitude  of  the  pressure  field  in  the  inner  fluid  as 


and 


O 

Piir,) 

—p—  =  • 

O 


(C-1) 


(C-2) 


The  circumferential,  longitudinal,  and  time  dependence  term  cos  (nO)  has  been 

dropped  from  the  equations  given  here.  The  values  calculated  are  at  a  maximum  amplitude; 
therefore  cos  (nO)  e  is  equal  to  1. 

The  following  equations  are  used  to  compute  displacements,  stresses,  and  strains  in  cylinder  2 
of  the  two-layer  cylinder.  If  the  cylinder  of  interest  was  the  first-layer  (rather  than  the  second 
layer  shown  here),  coefficients  \  \  cf  *  and  would  be  substituted  for 

C2  C2  C2  C2  C2  C2 

Aj  ,A2  ,5j  ,^2  and  C2  in  the  equations.  Quantities  from  the  single-layer  model  are 
obtained  by  using  coefficients  Aj  ,A2  ,  ^2  and  resulting  from  system  matrix 

(85). 
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The  magnitude  of  the  radial  displacement  of  a  particle  of  cylinder  material  is  given  by 
equation  (C-3)  as 


Cl 


^1  ^2  +  J-  (C-3) 


,C2. 


The  magnitude  of  the  longitudinal  displacement  of  a  particle  of  cylinder  material  is  given  by 


C2 


C2 

A]  iikJ^ip^fi))  +^2  (ikY^ip^r^))  + 


cf  (-1^.  ,(«2''l)  -  (^  K«  *  l<'?2''.)  J  + 


(C-4) 


The  radial  stress  in  the  cylinder  is  given  by 


C2  ,  , 

V  ("i) 


,C2 


A 


C2 


"  ^  X; 

(X,2  +  2^2)  “  '^l-^n^l^'l^ 

;)2  ^  ^ 

(^2  +  2II2)  ^^">2^1)  +  -  ^2^>2''l) 


-2  +  ^ 


y 


V'-i 

( 2 

■^  +  r 

yj 


A 


B 


Cl 


2[L^n 


-B 


Cl 


'2[i.nf  0 


+ 


Cj  2p2*^^‘^„+ l(^2^1^  ■*■  ^2  1  ^‘^2'’!^  J- 


(C-5) 
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The  magnitude  of  the  circumferential  strain  in  the  cylinder  is  given  by 


C2 


^00  (^i)  xif  13 


=  A 


B 


C2 


+ 


(C-6) 


Equations  (C-1)  through  (C-6)  have  been  normalized  by  P^.  If  a  different  excitation  was 

chosen  (i.e.,  or  Pq),  the  coefficients  \  \  \  \  cf  ^  and  cf  or 

C2  C2  C2  C2  C2  C2 

Ai  ,A2  ,5j  ,Cj  and  Cj  resulting  from  the  system  matrix  inversion  would  be  different, 
and  the  appropriate  excitation  would  be  substituted  for  in  equations  (C-1)  through  (C-6). 
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APPENDIX  D 

TWO-LAYER  MODEL  VERIFICATION 


In  this  section,  a  verification  of  the  two-layer  elasticity  model  relative  to  the  single-layer 
model  is  undertaken.  The  objective  is  to  produce  the  same  result  for  both  models.  First,  let  us 
consider  the  single-layer  cylinder  and  fluids,  whose  material  properties  and  dimensions  are  given 
in  tables  2,  4,  and  5  of  the  main  text.  The  evaluation  then  amounts  to  the  use  of  the  two-layer 
model  to  analyze  a  cylinder  that  is  physically  identical  to  the  single-layer  cylinder.  In  reality,  the 
two-layer  cylinder  is  homogeneous;  however,  it  is  treated  mathematically  as  a  two-layer 
composite.  The  material  properties  and  dimensions  of  the  two-layer  cylinder  used  for  the 
comparison  are  listed  in  table  D-1.  The  inner  and  outer  fluids  (tables  D-2  and  D-3)  have 
properties  consistent  with  those  of  tables  4  and  5. 


Table  D-1.  Two-Layer-Cylinder  Properties — 3.00-in.  Diameter  for  Model  Comparison 


Cylinder  #1  Property 

Cylinder  #2  Property 

Definition 

E  =  1.0x10^^ 

‘  m2 

II 

b 

X 

o 

00 

Young’s  Modulus 

Cl  =  0.3 

m 

b 

II 

Structural  Loss  Factor 

Pi  =  1070^ 

Pi  =  1070^ 

Density 

o 

II 

b 

II 

Poisson’s  Ratio 

a  =  1.200  in. 

b  =  1.350  in. 

Inner  Radius 

b  =  1.350  in. 

c  =  1.500  in. 

Outer  Radius 

Aj  =  0.15  in. 

=  0.15  in. 

Thickness 

D-1 
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Table  D-2.  Two-Layer-Cylinder  Outer  Fluid  Properties  (Water)  for  Model  Comparison 


Property 

Definition 

=  1000.0^ 

Density 

c  =  1500.0  ^ 

^  sec 

Speed  of  Sound 

Table  D-3.  Two-Layer-Cylinder  Inner  Fluid  Properties  (Water)  for  Model  Comparison 


Property 

Definition 

Pi  = 

Density 

C.  =  1500.0  ^ 

‘  sec 

Speed  of  Sound 

Simulations  are  performed  for  the  single-  and  two-layer  models  with  a  normal  pressure 
excitation  (n  =  0).  Figure  D-1  is  the  simulation  for  the  two-layer  model  and  figure  D-2  is  the 
simulation  for  the  single-layer  model.  In  figure  D-3,  a  comparison  of  the  two  models  is  made  at 
k  =  0  rad/m.  The  curves  in  this  figure  overlay  each  other,  which  indicates  agreement  between 
the  solutions.  Figure  D-1  is  subtracted  from  figure  D-2  and  the  result  is  displayed  in  figure  D-4, 
where  the  superscripts  M2,  and  Ml  correspond  to  the  two-layer  and  single-layer  models, 
respectively.  The  agreement  between  solutions  is  less  than  0.05  dB  everywhere  on  the  surface, 
except  at  very  low  frequencies  and  very  high  wavenumbers.  These  results  confirm  the  validity 
of  the  two-layer  cylinder  model  relative  to  the  single-layer  cylinder  model.  Figure  D-1  is  also 
used  on  the  cover  of  this  report. 

Figures  D-5  through  D-8  compare  the  two  models  (as  previously  described)  for 
circumferential  order  number  /z  =  1,  and  figures  D-9  through  D-12  compare  the  two  models  for 
circumferential  order  number  n  =  2.  In  both  cases,  the  agreement  between  the  models  is 
extremely  close;  no  measurable  difference  results  until  the  low-frequency/high-wavenumber 
region  is  reached. 
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re  D-1.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(ri)/P^)^  at  =  1.1  in. 
for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Two-Layer  Cylinder) 
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Figure  D-2.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P,(r|)/P^)^  at  =1.1  in. 
for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in.  (Single-Layer  Cylinder) 
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Figure  D-3.  Comparison  of  Two-Layer  and  Single-Layer  Elasticity  Models  With 
Magnitude  =  (10Log(P/ri;/Fo)^)  at  rj  =  1.1  in.,  Where  n  =  0  and  Diameter  =  3.00  in. 


Figure  D-4.  Relative  Difference  Between  Two-Layer  and  Single-Layer  Surfaces  With 
Magnitude  =  (10Log(P;^^(ri)/P„)2  -  at  =  1.1  in. 

for  Elasticity  Model,  Where  n  =  0  and  Diameter  =  3.00  in. 


m 


WAVENUMBER  (rad/m)  WAVENUMBER  (rad/m 


TR  11,067 


200.00 


100.00 


•100.00 


■200.00 


MAGNITUDE  (TWO-LAYER) 


20000.00  40000.00 

FREQUENCY  (Hz) 


60000.00 


re  D-5.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/ri)/P<,)^  at  =  1.1  in. 
for  Elasticity  Model,  Where  n  =  1  and  Diameter  =  3.00  in.  (Two-Layer  Cylinder) 
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Figure  D-6.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(P/r,)/P^)2  at  rj  =1.1  in. 
for  Elasticity  Model,  Where  n  =  1  and  Diameter  =  3.00  in.  (Single-Layer  Cylinder) 


1 


TR  11,067 


Figure  D-7.  Comparison  of  Two-Layer  and  Single-Layer  Elasticity  Models  With 
Magnitude  =  (10Log(P/r,)/Po)^)  at  ri  =  1.1  in.,  Where  n  =  1  and  Diameter  =  3.00  in. 


Figure  D-8.  Relative  Difference  Between  Two-Layer  and  Single-Layer  Surfaces  With 
Magnitude  =  (lOLogiPi^^irO/Po?  -  10Log(P;^Vri)/Fo)")  at  =  1.1  in. 
for  Elasticity  Model,  Where  n  =  1  and  Diameter  =  3.00  in. 
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re  D-9.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(F,<ri)/F„)2  at  rj  =  1.1  in. 
for  Elasticity  Model,  Where  n  =  2  and  Diameter  =  3.00  in.  (Two-Layer  Cylinder) 
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Figure  D-10.  Pressure  Transfer  Surface  With  Magnitude  =  10Log(Pj(rD/P<,)^  at  =1.1  in. 
for  Elasticity  Model,  Where  n  =  2  and  Diameter  =  3.00  in.  (Single-Layer  Cylinder) 
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Figure  D-11.  Comparison  of  Two-Layer  and  Single-Layer  Elasticity  Models  With 
Magnitude  =  (lOLogCP/rij/P^)^)  at  =  1.1  in.,  Where  «  =  2  and  Diameter  =  3.00  in. 


Figure  D-12.  Relative  Difference  Between  Two-Layer  and  Single-Layer  Surfaces  With 
Magnitude  =  {\m.og{Pi^\r^)IPo?  -  10Log(Pj^Vri)/PJ^)  at  =  1.1  in. 
for  Elasticity  Model,  Where  n  =  2  and  Diameter  =  3.00  in. 
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